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Abstract 

In this paper we study short time asymptotics of a density function of the solution of 



a stochastic differential equation driven by fractional Brownian motion with Hurst pa- 
rameter H G (1/2, 1) when the coefficient vector fields satisfy an ellipticity condition at 
the starting point. We prove both on-diagonal and off-diagonal asymptotics under mild 
additional assumptions. Our main tool is Malliavin calculus, in particular, Watanabe's 
theory of generalized Wiener functionals. 

Q . 1 Introduction 
(N 

Let {wt)t>o be the standard c?- dimensional Brownian motion and let (0 < i < d) be 
smooth vector fields on R" with sufficient regularity. Consider the following stochastic 
differential equation (SDE) of Stratonovich-type; 

X ■ dyt = ^ Vi{yt) o dw\ + Vfi{yt)dt with yo = a e R". 

; 

If the set of vector fields satisfies a hypoellipticity condition, the solution yt = yt{a) has a 
smooth density pt{a, a') with respect to Lebesgue measure on R". From an analytic point 
of view, pt{a,a') is a fundamental solution of the parabolic equation du/dt = Lu, where 
L = Vq + (1/2) Yli=i '^i 1 s-^d is also called a heat kernel of L. 

In many fields of mathematics such as probability, analysis, mathematical physics, 
and differential geometry, short time asymptotic of pt(a, a') is a very important problem 
and has been studied extensively. Although analytic methods are also well-known, we 
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only discuss a probabilistic approach via Feynman-Kac formula in this paper. Malliavin 
calculus is a very powerful theory and was used in many papers on this problem. 

Among them, S. Watanabe's result seems to be one of the best. (See [20] or Sec- 
tions 5.8-5.10, [8].) First, he constructed a theory of generalized Wiener functionals (i.e., 
Watanabe distributions) in Malliavin calculus. Then, he gave a representation of the heat 
kernel by using the pullback of Dirac's delta function; pt{a,a') = K[6a'{yt{a))], where the 
right hand side is the generalized expectation with respect to Wiener measure. Finally, 
by establishing an asymptotic expansion theory in the spaces of generalized Wiener func- 
tionals, he obtained a short time expansion of pt{a,a') under very mild assumptions. In 
this method, an asymptotic expansion is actually obtained before taking the generalized 
expectation. 

In this paper we consider the following problem. Let {w^)t>o be (i-dimensional frac- 
tional Brownian motion (fBm) with Hurst parameter H G (1/2, 1). Instead of the above 
SDK, we consider 

d 

dyt = ^ Vi{yt)dw^'' + Vo{yt)dt with yo = aeR"-. 

i=l 

This is an ordinary differential equation (ODE) in the sense of Young integral (see Lyons 
[12]). In fact, this is actually an ODE with a random driving path, but we call this 
SDE for simplicity. Some researchers have studied the solution of the above SDE with 
Malliavin calculus. See [161 E 113 HI E] and references therein. Under the ellipticity or 
the hypoellipticity conditon, the solution yt = yt{ci) has a smooth density pt{a,a'). See 

In this paper, by using Malliavin calculus and, in particular, Watanabe distribution 
theory, we will prove a short asymptotic expansion of this density in the elliptic case under 
mild assumptions. This kind of asymptotics was already studied in [H |2] , but without 
Malliavin calculus. In [1], they showed on-diagonal short time asymptotics when Vq = 0. 
In [2], by using Laplace's method, they showed off-diagonal short time asymptotics when 
Vq = and the vector fields V^'s satisfy a rahter special condition. Our results is a 
generalization of these preceding ones. Notice that we do not assume the drift term Vq is 
zero. One may think this is just a minor generalization, but this makes the asymptotic 
expansion much more complicated. 

The organization of this paper is as follows: In Section 2, we give settings, assumptions, 
and precise statements of two main theorems. In Section 3, we recall basic properties of 
a Young ODE and its Jacobian process for later use. In Section 4, we review Watanabe's 
theory of generalized Wiener functionals in Malliavin calculus. In Section 5, we discuss 
the solution of Young ODE driven by fBm with Hurst parameter H G (1/2,1) from 
the viewpoint of Malliavin calculus. We also prove uniform non-degeneracy of Malliavin 
covariance matrix of the solution under the ellipticity condition. In Section 6, we prove 
one of our main theorems, namely, on-diagonal asymptotics of the kernel. In section 7, we 
show the shifted solution of the Young SDE admits an asymptotic expansion in the sense 
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of Watanabe distribution theory. In Section 8, we prove the other of our main theorems, 
namely, off-diagonal asymptotics of the kernel. In Appendix, we make sure that Baudoin 
and Ouyang's result in [2] is basically included in ours. 

2 Setting and main results 

2.1 Setting 

In this subsection, we introduce a stochastic process that will play a main role in this pa- 
per. From now on, dropping the superscript "if", we denote by {wt)t>o = {w}, . . . , Wf)t>o 
the c?- dimensional fractional Brownian motion (fBm) with Hurst parameter H (1/2 < 
H < 1). It is a unique c?- dimensional mean- zero Gaussian process with covariance 

E[wlwi] = + -\t- sp^), (s, t > 0). 

Note that, for any c > 0, {wct)t>o and {c^Wt)t>o have the same law. This property is 
called self- similarity or scale invariance. 

Let Vi : R" — )■ R" be C^, that is, Vi is a bounded smooth function with bounded 
derivatives of all order {0 < i < d). We consider the following stochastic ODE in the 
sense of Young; 

d 

dyt = Y,Vi{yt)dwi + Vo{yt)dt with yo = a G R". (2.1) 

i=l 

We will sometimes write yt = ytia) = yt{ci,w) etc. to make explicit the dependence on a 
and w. 

2.2 Assumptions 

In this subsection we introduce assumptions of the main theorems. First, we assume the 
ellipticity of the coefficient of (12. ip at the starting point a e R". 

(Al): The set of vectors {Vi(a), . . . , Ki(a)} linearly spans R". 

It is known that, under Assumption (Al), the law of the solution yt has a density pt{a, a') 
with respect to the Lebesgue measure on R" for any t > (see [H [H]). Hence, for any 
measurable set U C R", P(?/t E U) = fjjPt{a,a')da'. 

Let H = be the Cameron-Martin space of fBm (wt). For 7 G we denote by 
0° = 0f (7) be the solution of the following Young ODE; 

d 

d(j)'i = J2Vi{(pt)d-fi with 0° = aGR". (2.2) 

1=1 
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Set, for a 7^ a', 

< = {ieH I 0°(7) = a'}. 

If we assume (Al) for all a, this set K^' is not empty. If K^' is not empty, it is a Hilbert 
submanifold of Ti. From the Schilder-type large deviation theory, it is easy to see that 
inf{||7||^ I 7 G K^'} = min{||7||-^ | 7 G K^'}- Now we introduce the following assumption; 

(A2): J eK^ which minimizes H-norm exists uniquely. 

In the sequel, 7 denotes the minimizer in Assumption (A2). We also assume that || • ||^/2 
is not so degenerate at 7 in the following sense. 

(A3): At 7, the Hessian of the functional K^' 9 7 i— )■ ||7|||^/2 is strictly larger than 
Id^ff/2 in the form sense. More precisely. If (— ^o^^^o) 3 m >— ?■ f{u) E K^' is a smooth 
curve in such that /(O) = 7 and /'(O) ^ 0, then {d/du)\=o\\f{u)fn/2 > ||/'(0)||2^/2. 

Later we will give a more analytical condition (A3)', which is equivalent to (A3) under 
(A2). In [20], Watanabe used (A3)'. We will also use (A3)' in the proof. In order to 
state (A3)', however, we have to introduce a lot of notations. So, we presented (A3) 
here for convenience. 



2.3 Index sets 

In this subsection we introduce several index sets for the exponent of the small parameter 
e > 0, which will be used in the asymptotic expansion. Unlike in the preceding papers, 
index sets in this paper are not the set of natural numbers and are rather complicated. 
Set 

Al = {ni + I ni,n2 E N}, 

where N = {0, 1, 2, . . .}. We denote by = < < K2 < ■ ■ ■ all the elements of Ai in 
increasing order. Several smallest elements are explicitly given as follows; 

1 12 

Kl = 1, '^2 = Jj, % = 2, K4 = l + — , K5 = 3A— 

As usual, using the scale invariance (i.e., self-similarity) of fBm, we will consider the scaled 
version of (12.11) . (See the scaled Young ODE (16.11) below). From its explicit form, one can 
easily see why Ai appears. 
We also set 

A2 = {«:-1 I «:eAi\{0}} = {0, -i-1, 1, ^, (3 A |:) - 1, . . .} 

and 

A', = {k-2\keAi\ {0, 1, 1/H}} = jo, 1 - 1, (3 A -|) - 2, . . .}. 
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Next we set 



A3 = {ai + a2 + ■ ■ ■ + a^i \ m ^ N+ and ai, . . . , G A2}. 

In the sequel, {0 = z/q < ^1 < ^2 < ■ ■ ■ } stands for all the elements of A3 in increasing 
order. Similarly, 

A3 = {ai + a2 + ■ ■ ■ + ttm \ m E N+ and ai, . . . , G A'2}. 

In the sequel, {0 = po < Pi < P2 < ■ • ■ } stands for all the elements of A3 in increasing 
order. Finally, 

A4 = A3 + A3 = {z/ + p I z/ G A3, p G A3}. 
We denote by {0 = Aq < Ai < A2 < ■ ■ • } all the elements of A4 in increasing order. 

2.4 Statement of the main results 

In this subsection we state two main results of ours, which are basically analogous to the 
corresponding ones in Watanabe |20]. However, there are some differences. First, the 
exponents of t are not (a constant multiple of) natural numbers. Second, cancellation of 
"odd terms" as in p. 20 and p. 34, [20j does not happen in general in our case. (If the 
drift term in Young ODE (12. ip is zero, then this kind of cancellation takes place as in 

mm)- 

The following is a short time asymptotic expansion of the diagonal of the kernel 
function. This is much easier than the off-diagonal case. 

Theorem 2.1 Assume (AI). Then, the diagonal of the kernel p{t,a,a) admits the fol- 
lowing asymptotics as t \0; 

p{t, a,a)^^ (co + c,,t"i^ + c,,r^^ + ■ ■ ■ ) 

for certain real constants Cq, c^^, Cj,2, • • •• Here, {0 = z/q < ^1 < ^2 < ■ ■ ■ } (if^G all the 
elements of A3 in increasing order. 

We also have off-diagonal short time asymptotics of the kernel function. 

Theorem 2.2 Assume (A1)-(A3). Then, we have the following asymptotic expansion 
as t \ 0; 

Pit, a, a') ~ exp + ^) ^{"Ao + «A,t^^'' + ax.t^''' + ■■■} 

for certain real constants (3, (j = 0, 1, 2, . . .). Here, {0 = Aq < Ai < A2 < ■ ■ • } are all 
the elements 0/A4 in increasing order. 
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Remark 2.3 (i) Consider the following simplest case; n = d = 1 and yt = a + wt + bt 
with 6 G R. Then, for each t > 0, this induces a Gaussian measure with mean a + bt and 
variance t^^ . Hence, the kernel is given by 



pit, a, a') = , — exp ( 



1 / (a + 6t — a' 



/\2 



1 „l\2llot2H\ _h(„_„'\/^2H-l u2,2-2H 



^-{a-a')y{2t^")^-b{a-a')/t^H-i^_^2t2-2H/^ 



2'Kt" 

{a-a'f /{2t^H)-b{a-a')/t^"-^ , 



^t^\ 2' +222!' )■ 



This example may illustrates that asymptotics in Theorem is not so strange. 

(ii) Some of the constants in Theorems \2.1\ and \2.2\ can be obtained explicitly. For 

example, in Theorems \2.1\ Cq = [(27?)"/^ det(a(a)cr(a)*)]~-'^ and 

n 

c^i = C(i/H)-i = + ■ ■ ■ + Ma)wt) ■ Vo{ay = 0. 

Here, o'{a)a{a)* is the covariance matrix of the n-dimensional Gaussian random variable 
Yl'j=i^ji^)'^i- -^'^ Theorems \2.2[ (5 = {u^cjy^^). Notations in this remark will be given 
later. 



Basic properties of Young ODE and L^-integrability 
of Jacobian process 



In this section we recall the basic properties of a Young ODE and its Jacobian process 
(i.e., derivative process). These facts scatter in many literatures and it is not so easy to 
find a suitable one. (In this sense, Lejay [10] may be useful.) Here, we summarize some 
results for later use. 

Most of the contents in this section are not new. The only new part is L'^-integrability 
of the Jacobian process driven by fBm with Hurst parameter H > 1/2 for any q > 1. 
(This may also be known to the experts, however, although there seems to be no literature 
that explicitly proves it.) 

We always assume that 1 < j9 < 2 and the time interval is [0, 1]. Let C^^p-''^'^{[0, 1]; R'^) 
and C^~'""^'([0, 1]; R*^) be the spaces of R"'- valued Holder continuous paths and con- 
tinuous paths of finite p-variations, respectively. The Banach norms are defined by 



\x\\i/p-hid = \xo\ + sup 



o<s<t<i {t-sy/p' 

l^^lll/p-mr = l^^ol + sup ( y \Xt^ — Xt^_^\^ I 



1 



i=l 
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where sup runs over all the finite partition P = {0 = to < < ■ ■ ■ < ^iv = 1} of [0, 1]. 
The closed subspaces of paths that starts at the origin is denoted by Co/^"'''''([0, 1]; R"^) 
and ^^-''"'■([0, 1]; R'^), respectively. 

Let a : R" — )■ Mat(n, c?) be sufficiently regular. Consider the following ODE in the 
Young sense; 

dyt = a{yt)dxt with = a. (3.1) 

Here, x is a continuous path in R"^ of finite p-variation and a G R" is the initial value. 
Let : R" R" be the ith column vector oi a (1 < i < d). Then, ODE (EH]) can be 
rewritten equivalently as follows; 

d 

dyt = ^Viiytjdxl with yo = a. (3.2) 

1=1 

Some researchers prefer this style. In this paper we will use both (13. ip and (13. 2p . 

If a is C|, then the above ODE has a unique solution for any given x and a in p- 
variational setting. Moreover, the map 

C^"'^'^([0, 1]; R'^) X R" 9 (x, a) H-> y G C"^'([0, 1]; R") (3.3) 

is locally Lipschitz continuous (i.e., Lipschitz continuous on any bounded set). 

To be more precise, the following well-known propositions hold. (It is important that 
the estimates in Propositions 13. II and 13.21 are independent of the initial value a.) In what 
follows, p G [1,2) is fixed. We set K = maxo<i<2 ||VV||oo, where || ■ ||oo denotes the 
sup-norm. 

Proposition 3.1 Let a be Cl and x G C^^™'^([0, 1]; R'^) and consider ODE (EJ]) in the 
Young sense. Suppose that a control function uj satisfy that 

\xt -Xs\< w(s, tf'P (0 < s < t < 1). 

Then, the following (i), (ii) hold: 

(i) There exists Ti G (0, 1] such that, on the restricted interval [0,Ti], a unique solution 
y of 113. 1\) exists and satisfies the following estimate; 

\yt-ys\<2Ku{s,ty/P (0<s<t<Ti). 

(Ti can be chosen as follows. For a certain constant C G (0, 1], which is independent of 
the initial value a, we may set Ti = sup{t G (0, 1] | uj{0,tY^'P < C/(l + K)"^}. Then, the 
above estimate holds.) 

(ii) In the same way as in (i), we have a partition = Tq < Ti < ■ ■ ■ < T^ = I of 
[0, 1] and a unique local solution on each subinterval [Tj_i,Tj]. Concatenating these local 
solution, we obtain a global unique solution y, which satisfies the following estimate; 

b. - ,,| < 2k{i + (0 < . < t < 1). 

where R = C/{1 + K)\ 
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Proposition 3.2 Let a be Cl and x,x e ^^"^"''([0, 1]; R'^) and consider ODE / TO) in 
the Young sense. Suppose that x,x, e > 0, and a control function u satisfy that 

\xt — Xs\, \xt — Xs\ < uj{s,ty^^ (0 < s < t < 1), 

\{xt-Xs)-ixt-Xs)\ <eu{s,ty/P (0<s<t<l). 

Then, the following estimate holds; 

\{yt - ys) - [m - ys)\ < e[l + ^)^'"^''2^+-(°'^)/«'a;(.,t)V. (0 < . < t < 1). 

Here, R! = C'/{1 + K)'^ for a certain constant C G (0,C], which is independent of the 
initial value a. 

Remark 3.3 Propositions \3.1\ and \3.^ cover the Holder case, too. For example, if x is 
l/p-Holder continuous, we may takeu{s,t) = W^Wi/p^hifii't — sY^^ in Proposition \3 . 1\ (i). 
In particular, local Lipschitz continuity of the ltd map as in Ii3.3\) holds for l/p-Holder 
continuous case, too. 

Now we discuss the Jacobian process (i.e., the derivative process) J of the ODE (13.11) . 
or equivalently (13.21) . Jt is a (formal) derivative of the solution flow yt = yt{ci) of the 
Young ODE (l33|) . 

For V G R", we denote the directional derivative along v by V^,c^(J/) = Vo"(?/)(w, ■ ), 
etc. So, V(T takes its values in L(2)(R", R'^; R'^) = (R")* ® (R'^)* ® R", which is equipped 
with the usual Hilbert-Schmidt norm. Notations such as V*l^-, V^cr = VVcr, etc. should 
be understood in a similar way. 

The Jacobian process J takes its values in Mat(r2, n) = L(R", R") and satisfies 

dJt = V(T{yt){Jt,dxt) with Jq = Id„. (3.4) 

More precisely, by setting Mt = /J Vcr(ys) ( ■ , dxg) , we may rewrite this equation as follows; 

dJt = dMt ■ Jt with Jo = Id„. (3.5) 

The dot on the right hand denotes the matrix multiplication. When we need to specify 
the driving path, we will write J{x). The equivalent equation for J that corresponds to 
(13. 2p is as follows; 

d 

dJt = J2^V,{yt){Jt)dxl with Jo = Id^„. (3.6) 

i=l 

It is known that the system of Young ODEs (13. ip and (13. 4p has a unique solution 
{y, J). We need a nice estimate for J in order to prove that Holder norm of J has all 
finite moments when the driving path x is fBm with H > 1/2. Hu and Nualart [7] proved 
that supo<«i I Jt| G nq>iL^, but we will estimate Holder norm of J below. 
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Since (13 ■4p is linear, the solution can be written explicitly as follows. 

oo 

Ji= (ld„ + ^M]5)j, (0<s<t<l), (3.7) 



A;=l 

where 



Mi'^l = / dMt^---dMt,dMt,. (3.8) 

J0<tl<■■■<t^:<l 

By a special case of Lyons's extension theorem (p. 35, [13]), we have the following 
estimate (13.91) below. Assume that M satisfies |m]\'| < a;(s, t)^/^//3(l/p)! for a control 
function w on a certain times interval [a,b], where /3 > is a constant such that > 

l + Etsilhf''"''^''- Then, we have 

'^''^'-T^MT' {a<s<t<b, keN^). (3.9) 

Let Ti be as in Proposition 13.11 On this interval, 

\xt -x,\ + \yt-y,\< {2K + l)u{s, tf'^ (0 < s < t < Ti). 
Moreover, the following map from R" to L(R'^, Mat(n, n)) 

R" 3 y ^^ [x e R"' ^-> V(T(y)(-,x) G Mat(n,n)] 

is bounded and globally Lipschitz continuous with Lipschitz constant K. Hence, by 
the theory of Young integral, there exists a positive constant ci^k = ci,K,p such that 
|m]5| = \Mt - Ms\ < Ci,i^w(s,t)i/P for all < s < t < Ti. If we choose such a /3 > 
and set a)(s,t)^/^ = ci,ii:/9(l/p)! ■ u{s,ty^P, we may use Lyons's extension theorem for 
this control function u and obtain (13. 9p . Note that a)(0,Ti)^''^ < C2,k for some constant 
C2,K = C2,K,p > which may depend only on K, p. 
It is easy to see that 



El /Iff'"' I < '^y^^'^) ^ ujyt,,i) >^ 



k=l k= 



Consequently, supo<t<Ti \ Jt\ < 1 + C3,_ft:C/(l + i^T)^ = c^^k and 

oo 

l^i - ^.1 < l^sl J2 I^Sl ^ C4,i^C3,xC^(s,t)^/P, (0 < S < t < Ti). 

k=l 

Here, Cs^x and C4^x may depend on K,p, (and the choice of (3), but not on the initial value 
a. 
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Let = Tq < Ti < T2 < ■ ■ ■ < Tn = 1 he a partition of [0, 1] such that 
w(0, Ti) = ■ ■ ■ = u(Tn-2, Tn-i) = R> uj(Tn--i, 1), 

where R = C/{1 + KY. Note that — 1 < u;(0, holds since u is super additive. 
Then, on each subinterval [Tj_i,Tj], we have Proposition 13. (i) and hence, we have the 
inequahty f l3.10p for Tj_i < s < t < Tj. Note that we can take the same constant c^^k- 
By induction, we can easily see that sup2-^_^<j<j.^, \Jt\ < c^/^ and that 

00 

\Jt - Js\ < \Js\ J2 < C^,,KC3,Ms,ty/P, (r,_i <S<t< T,). 

k=l 

Therefore, on the whole interval, we have 

\Jt - Js\ < CW(1 + ^^Y'~'"Ms.tfl^ 

< C3,^exp[(l + ^)logC4,^] (1 + ^)(-^)/-c.(.,t)V. (3.11) 

for any < s < t < 1. This shows that p- variation norm of J is of exponential order of 
c<;(0, 1). In particular, if x is 1/p-Holder continuous, we can take uj{s,t) = \\x\\yp_f^i^(t — s) 
and we have from (13. lip that 

\\J\\i/p-hid < 1 + C5,k{1 + Ml/p^hid) exp(c6,i^||x||^/p_^J. (3.12) 

Here, positive constants c^^k and cq^k depend only on K,p, (and the choice of (3). 

Next we consider J^^. Instead of (13. Sp . (13. 7p . and (13. Sp . we have similar equations as 
follows; 

dJ-^ = -J^^-dMt with J(7^ = Id„, (3.13) 

00 

J-i = J7'(ldn + $^M]?), (0<s<t<l), (3.14) 

k=l 

where M^} = {-iff dMt.dMt, ■ ■ ■ dMt^. (3.15) 

From these, we can easily see that the estimates (13. lip and (I3.12p holds for J~^, too. 

Now let us consider the moment estimate for J and J~^, when the driving path x 
is the (i-dimensional fBm w = (wt)o<t<i with Hurst parameter H G (1/2, 1). Take any 
p G (1/if, 2). Then, almost surely, Hwlli/p-ziid < 00. (By the way, = 00, a.s. 

See [S1[IB].) 

The differential equations are given as follows; 
dyt = a{yt)dwt with yo = a and dJt = dMt ■ Jt with Jq = Id„, (3.16) 

where Mt = Jl^Vcr(ys){- ,dws)- For simplicity we call it an SDK, though it is just a 
deterministic Young ODE driven by a random input w. 
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Proposition 3.4 Let H and p be as above and assume that the coefficients a is C^. Let 
J be as in Ii3.16\) above. Then, \\J\\i/p-hid one? have moments of all order, 

i-e., \\J\\i/p-hid, \\J^'^\\i/p-hid G ni<g<oo-^^'^- 



Proof. This is immediate from (13 .121) and Fernique's integrability tlieorem for Gaussian 
measures. I 



Remark 3.5 It is also possible to prove Proposition \3.4\ by using Hu-Nualrt's moment 
estimate for suPq<^<;^ | Jj| in ^ and a cutoff argument. 

Remark 3.6 This is a remark on a Young ODE with a drift term. Let b = Vq : R" — )■ R"" 
be Cl. Instead of h3. 1\) or hS. we consider 

dyt = a{yt)dxt + h{yt)dt with yo = a, 

or equivalently, 

d 

dyt = ^ Vi{yt)dx\ + Vo{yt)dt with yo = a. 

i=l 

All the results in this section hold true with a trivial modification for a Young ODE with 
a drift term above. The proof is easy. We just need to replace d,x,a with d + l,x,a, 
respectively. Here, x is a {d + 1)- dimensional "space-time path" defined by xt = {xt,t) 
and a is an n X (d + 1) -matrix-valued function defined by [a\b] = [Vi, . . . , V^, Vq]. 
In this case, and l{3.6\) read 

dJt = Va{yt){Jt,dxt) +Vb{yt){Jt)dt with Jo = Idn, (3.17) 

d 

dJt = ^Vi{yt){Jt)dxi + VVo{yt){Jt)dt with Jo = Idn, 

i=l 

respectively. When necessary, we will write J = J{x, X), where Xt = t. 

4 Preliminaries from Watanabe's asymptotic theory 
of generalized Wiener functionals 

In this section we recall Watanabe's theory of generalized Wiener functionals in Malliavin 
calculus. Most of the contents in this section are borrowed from Sections 5.8-5.10, Ikeda 
and Watanabe p]. Shigekawa [19] and Nualart [15] are also good textbooks of Malliavin 
calculus and we will sometimes refer to them. There is no new result in this section. 
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First, let us recall the basic notions in Malliavin calculus. Let {W, "H, /i) be an abstract 
Wiener space. As usual, "H and "H* are identified through Riesz isometry. Let D be the 
'H-derivative and D* be its dual. Ornstein-Uhlenbeck operator is denoted by L = —D*D. 
The first Wiener chaos associated with h E Ti is denoted by {h,w). If {h, ■) E W*, then 
it coincides with ]v*{h,w)w A function F : — )■ R is called a real-valued polynomial if 
there exist m G N, hi, ... , hm, and a real- valued polynomial / of m- variables such that 
F{w) = f{{hi,w),...,{hm,w)). The set of all real-valued polynomials are denoted by 
V. Let /C be a real separable Hilbert space. U{1C) = L'^(W;IC) denotes the L'^-space of 
/C-valued functions. A function G : W ^ K, is called a /C-valued polynomial if there exist 
m G N+, Fi, . . . , Fjn G V, and Vi, . . . ,Vm G JC such that G{w) = X]j=i Fj{w)vj. The set 
of all /C-valued polynomials are denoted by V{JC). 

For q G (l,oo) and r G R, and F G P(/C), set = ||(/ - L)-''^^F\\g. We define 

the Sobolev space Dg ,,(/C) to be the completion of V{IC) with respect to this norm. When 
/C = R, we simply write L'^, Dg^^, etc. If g < g' and r < r', then || ■ \\q^r < || ■ \\qiy and 

When r = G N, then the well-known Meyer equivalence holds (see [19] for example); 
there exists a positive constant Cg^r such that 

c;i\\F\\g,k < \\F\\g + WD'^Fllg < Cg,k\\F\\g,k for all F G Dg,fc(/C). 

Now we recall generalized Wiener functionals, which are often called Watanabe dis- 
tributions, too. Set Doo(/C) = ni<g<oo,rGRDg,r(/C) and D_oo(/C) = Ui<g<oo,rGRDg,r(/C). 

Those are called the space of test functions and the space of generalized Wiener function- 
als, respectively. We also use the following spaces of (generalized) Wiener functionals: 
Doo(A:) = n^^i Ui<g<oo D,,fc(/C) and D_oo(/C) = U^^i ni<g<oo D,,_fc(/C). In other words, 
F G Doo(/C) is equivalent to that, for any A; > 0, there exists q = q{k) > 1 such that 
F G D,,fc(/C). 

Let F = (F^,...,F") G Doo(R")- The Malliavin covariance matrix is defined by 
{^{DF^iw), DF^iw))^) We say that F is non-degenerate in the sense of Malliavin 

if det((DF*(w), D F\w)) T-i) ^l^. .^^ is in L'^ for all 1 < g < oo. This non-degeneracy is very 
important in Malliavin calculus. For example, such a non-degenerate F can be composed 
with a Schwartz distribution (p defined on R" and (p o F becomes a generalized Wiener 
functional. 

We introduce seminorms on the space of tempered distributions on R". For A; G Z and 
a real- valued, rapidly decreasing, smooth function (f) of Schwartz class 5(R") on R", we 
define \\(j)\\2k = ||(1 + I ■ P - A/2)V||oo. Set 52fc(R") to be the completion of 5(R") with 
respect to this norm. Then iS(R") = n,fc>oiS2A:(R"), which is a Frechet space. The dual 
space iS'(R") = U,fc>oiS_2fc(R") is called the Schwartz space of tempered distributions. 

S. Watanabe proved that, if a Wiener functional F is non-degenerate, then the pullback 
map I— >■ o F extends to a continuous linear map between two distribution spaces. The 
following is borrowed from pp. 378-379, Ikeda and Watanabe [8]. 
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Theorem 4.1 Let F = {F^, . . . , F") G Doo(R") be non- degenerate in the sense of Malli- 
avin. Then, for any 1 < g < oo and k = 0,1,2,..., there exists a positive constant 
C = C{q, k, F) such that 

||0oF||,,_2, <C||0||_2A. (0G5(R")) 

holds. Therefore, the map i— )■ o F extends uniquely to a continuous linear map 
5_2fe(R") BT^ToFe D,,_2fc. In particular, ToFe U^^^ ni<g<oo Dg,_2fc = D_«,. 

We recall an integration by parts formula in the context of Malliavin calculus (See p. 
377, [8j). For F = {F\ . . . , F") G Doo(R"), we denote by t'^{w) = {DF%w), DF^{w))n 
the (i, j)-component of Malliavin covariance matrix. We denote by 7*-'(w) the {i,j)- 
component of the inverse matrix r"^. Note that r*-' G Dqo and D^-' = / 7*'^(-D'r'^')7'-'. 
Hence, derivatives of 7*-^ can be written in terms of 7*-''s and the derivatives of r^-^'s. 
Suppose G G Doo and T G 5'(R"). Then, the following integration by parts holds; 

E[{d,T oF)-G]=E[{ToF)- $,,( ■ ; G)] (4.1) 

where $j( ■ ; G) G Doo is given by 

d d 

^i{w-G) = E G{w)^'\w)^^\w){Dt'\w),DF^{w))u 

3=1 k,l=l 

+ Y\w){DG{w),DF^{w))n + j'^{w)G{w)LF^{w)y (4.2) 

Note that the expectations in (14. ip are the pairing of D_oo and Doo- 

Now we consider a family of Wiener functionals indexed by a small parameter e G (0, 1]. 
Consider a family of /C-valued Wiener functionals {F(£, w)}o<£<i and assume F{e, ■ ) G 
Doo(/C) for each e. We say F{e, ■ ) = 0{e'^) in Dg_A:(/C), k G R, as £ \ 0, if ||F(e, ■ )\\g^k = 
Oie"). We say F{e, ■ ) = 0{e^) in Doo(/C) as e \ 0, if F{e, ■ ) = 0(8") in Dp,fc(/C)' for 
any 1 < g < 00 and A; G N. 

Let = Ko < < ft;2 < ■ ■ • 00 and /o, /^j, /«;2, • • • G Doo(^)- We write 

F(£, 0-/0 + e"'f.^ + + • • • in Doo(/C) as e \ 0, 

if, for any m G N, it holds that 

F{e, ■) - (/o + /., + ■■■ + e^'-f.J = 0(5'^'"+^) in Doo(/C) as £ \ 0. 

In a similar way, we define the asymptotic expansion in the space of Watanabe distri- 
butions. Suppose that $(£, ■ ) G D_oo(A^) for all < e < 1 and 0o, 0ki, 0ft2! • • • G D_oo(A^)- 
We write 

$(£, ■ ) ~ 00 + ^'''^^i + £'^'0.2 + ■ ■ ■ in D_oo(/C) as £ \ 0, 
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if, for any m = 0, 1, 2, . . ., there exist G N and 1 < g < oo such that 

■ ), 0,^ G D,,_fc(/C) (0 < £ < 1, < J < m) 

and 

■ ) - (00 + ^"V.i + • ■ • + = 0(5'^-+^) in D,,„,(/C) as £ \ 0. 

Let G{e, ■ ) G Doo(A^) for all < e < 1 and go, g^^,g^^,... G Doo(^)- We write 
■ ) ~ ^0 + e^^g^^ + e'^^c/^^ + • • • in Doo(/C) as e \ 0, 
if, for any m G N and /c G N, there exists 1 < g < oo such that 

■ ), g^^ G Dg,fc(/C) (0 < e < 1, < j < m) 

and 

■ ) - (^70 + + ■ ■ ■ + e'^-g^J = 0(5'^-+^) in D,,fc(/C) as e \ 0. 

Finally, let ^(e, ■ ) G D_oo(/C) for all < e < 1 and ?/'o, V'ki, ^^2? • • • ^ D-oo(/C). We 
write 

^!{e, ■ ) ~ V'o + ^'''V'^i + e'^'^Kz + • ■ ■ in D-oo(/C) as e \ 0, 
if, for any m G N, there exists A; G N such that 

^'(e, ■ ), e ni<g<ooDg,_fc(/C) (0 < £ < 1, < J < m) 

and, for any 1 < q < oo, 

^{e, ■ ) - (V^o + ^"^^.1 + ■ ■ ■ + e^-ij^J = Oie^-+') in Dg,_fc(/C) as e \ 0. 

We recall basic facts for such asymptotic expansions in the Sobolev spaces. Let = 
kq < ft^i < < ■ ■ ■ C)0 be as above. In Proposition 14. 21 below. = < ui < 1^2 < ■ ■ ■ 
00 are all the elements of {ni + kj | i, j G N} in increasing order. The fundamental case 
Kj = j is treated in Proposition 9.3, Section 5.9, Ikeda and Watanabe [H]. The following 
is a straight forward modification of it. 

Proposition 4.2 (i) Suppose that F{e, ■ ) G Doo(/C) admits an expansion such as 

F{e, ■ ) ~ /o + e^'f^, + e'^'f^, + ■■■ m B^ilC) as e\ 0, 

with £ Doo(A^) for all j G N. Suppose also that G{e, ■ ) G Doo (or Dooy' admits an 
expansion such as 

G{e, ■ ) ~ fi'o + ^""^Sk^ + ^"^9^2 H Doo (or resp. t)oo) as e\0, 
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with g^. G Doo (or resp. Tl>oo) for j E N. Then, H{e,w) = F{e,w)G{e,w) satisfies that 

H{e, ■ ) ~ ^0 + e^^Ki + £"^^2 H Doo(/C) (or resp. Doo(/C) ) as e\0, 

where hy^ G Doo(/C) (or resp. ti^iJC)) are given by the following formal multiplication; 

^l^n ~ ^ ^ dnifKj- 

{i,j);Ki + Kj=Un 

(ii) Suppose that G{e, ■ ) G Dqo (or T)oo) admits an expansion such as 

G{e, ■ ) ~ fi'o + e'^^gni + e^'^'fl'^a H in (or resp. Boo) as e\0, 

with g^^ G Doo (or resp. Dooj for all j G N. Suppose also that • ) G D_oo(/C) admits 
an expansion such as 

■ ) ~ 00 + + £"¥^2 + • • ■ D_oo(/C) as £ \ 0, 

wt/i G D_oo(A^) /o?^ fl^O ^ N). r/ien, = G{e,w)^{e,w) satisfies that 

^!{e, ■) r^i)Q + e'''i)u,+e''^i)u2 + --- mt>^^{lC) (or resp. D_oo(/C)) as e \ 0, (4.3) 

where ip^n ^ D„oo(A^) ^or resj>. D_oo(A^) y' o^^e g'Zi'en 6y the following formal multiplication; 

(iii) Suppose that G{e, ■ ) G Doo admits an expansion such as 

G{e, ■ ) ~ 5'o + + e'^^gK2 H Doo as e \ 0, 

with g^. G Doo for all j G N. Suppose also that ■ ) G D_oo(^) admits an expansion 
such as 

■ ) ~ 00 + ^'''^Ki + e''^<pK2 + • • • «^ D_oo(/C) as £ \ 0, 
with (pK. G D_oo(A^) for all j G N. T/ien, (^T^ anc? ( (^.^[ j /ioM m D_oo(A^)- 



Remark 4.3 /n (i) of the above Proposition, the index sets {kj}j=o,i,2,... for the asymp- 
totic expansions for F{e, ■ ) and G{e, ■ ) are the same. However, these index sets for F 
and G may differ, because the union of the two index sets can be regarded as a new index 
set. Similar remarks hold for (ii) and (iii), too. 
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Next we consider asymptotic expansions for the pullback. Let F{e, ■ ) G Doo(R"') for 
< £ < 1. We say F is uniformly non-degenerate in the sense of Malhavin if 

hmsup II det{{DF\e, ■ ),DF^{e, ■ ))h)7<, .<„||g < oo for all 1 < g < oo. 

The following is a straight forward modification of Theorem 9.4, [8]. In this theorem, 
= z/q < z/i < z/2 < • ■ ■ oo are all the elements of 

{/tji H V I n = 1, 2, . . . , and ji, . . . , j„ G N} 

in increasing order. 

Theorem 4.4 Let F{e, ■ ) G Doo(R."') (0 < £ < 1) satisfy the following; 

F{e, ■ ) ~ /o + e^'f., + s^'f.2 + " " " B^{R^) ase\0, 

with /kj. G Doo(R.") for all j G N. We also assume that F is uniformly non- degenerate 
in the sense of Malliavin. Then, for any T G iS'(R"), ^{e,w) := T o F{e,w) has the 
following asymptotic expansion; 

■ ) ~ 00 + e'"'(pK, + e^'^Ka + ■ ■ ■ in D_oo ase\0, 
where ^ D_oo is determined by a formal Taylor expansion as follows; 

■) = y2 ^id'^T)ifo)[s^'f., + e^^f^, + ...]- = 00 + e''^<P,, + e'^^0., + ■ ■ ■ , 
^ — ^ al 

a 

where the (formal) summation is over all multi-indexes a = [ai, . . . , a„) G N". (We set 
(9" = Ylji^/dx^)"^ and b" = Ylj^'j' /^'^ ^ — ipi^ ■ ■ ■ i^n) ^ R-" o^s usual.) For instance, 
00 = T(fQ) and 0«,, = YTj=i fL ' idT/dx^){fo) and so on. 

Unlike in the usual stochastic analysis, almost every Wiener functional in this paper is 
continuous with respect to the topology of an abstract Wiener space, because we work in 
the framework of Young integration. Therefore, the following proposition will be very use- 
ful. For Banach spaces Xi,...,Xm,y, L^"") (A'l, . . . , Xm] y) denotes the space of bounded 
m-multilinear maps from Xi x ■ ■ ■ x to y . 

Proposition 4.5 Let {Wj'Hjfi) be an abstract Wiener space. Then, we have the following 
bounded inclusions; 

L^"'\ W,..^,W ]R) ^ L^"'\ W,..^,W ,n;R) ^ [n*)^"". 

m m—1 

Here, the tensor product on the right hand side is Hilbert- Schmidt as usual. 

Proof. The left bounded inclusion is obvious. The right one is in p. 103, Kuo [9]. I 
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5 Some results on Malliavin calculus for the solution 
of Young ODE driven by fBm with H > 1/2 

In this section we discuss the solution of Young ODE driven by fBm with Hurst parameter 
H G (1/2,1). We give moment estimates for the derivatives of the solution and prove 
uniform non-degeneracy of Malliavin covariance matrix of the solution. 

Take a G {1/2, H) and set p = 1/a. We denote by /i = fi^ the law of (i-dimensional 
fBm starting at 0. This Gaussian measure is supported in ^^"'^''^([0, 1]; R'^). Cameron- 
Martin space is denoted hy n = , which is not dense in Cq-^^\[Q, 1]; R'^). We set W 
to be the closure of U in C^^''''^([0, 1]; R"'). Then, (>V, fi) becomes an abstract Wiener 
space. (Note that the separable Hilbert space n is not dense in C^-^^'^{[0, 1]; R'='), which 
is not separable.) We denote by {wt)o<t<i = {'wf)o<t<i the canonical realization of fBm. 

From now on, we assume that a : R" — > MaX{n,d) and b : R" — )■ R" are C^. We 
recall Young SDE (12. ip driven by fBm (wt) in the following form; 

dyt = a{yt)dwt + b{yt)dt with yo = a. (5.1) 

Then y{w) = I{w, A), where Xt = t and I is the Ito map corresponding to the coefficients 
[a; b] = [Vi,..., Vd, Vq]. I is everywhere-defined and continuous from Cg "''''^([0, 1]; R''^^) 
to C"^'^''^([0, 1]; R''), as we have explained in Section [31 

Moreover, I is smooth in Frechet sense (See Li and Lyons [llj) and, in particular, 
y = J( ■ , A) is infinitely differentiable in ?^-direction (see Nualart and Saussereau p!7]). 
These are deterministic results. In the sense of Malliavin calculus, it is shown in Hu and 
Nualart that yr : W ^ W is for any T e [0, 1]. 

We can obtain an explicit form of the directional derivative := Dhyt {h G "H) by 
differentiating (15. ip : 

rfef - Va(|/i)(ef , dwt) - Vb{yt){i'l)dt = a{yt)dht with = 0, (5.2) 

or equivalently, 

= J{w, \)t [ J{w, X);'a{yt)dht. (5.3) 

Note that all the integrations above are in the Young sense. An explicit form of J = 
J{w, A) is given in (13.171) . Let h,k E Ti. By differentiating the above ODE, we see that 
^k,h ._ jj^D^y^ satisfies the following ODE; 

d^t'' - ^<y{yt) {it'\ dwt) - vb{yt) {^^•')dt = V V(y,) (ef , , dwt) 

+ Va{yt){^t dkt) + Va{yt){^^, dht) + V%{yt){^t ^^)dt with ^o'" = 0. (5.4) 
Equivalently, we have 

^^/ = J{nj,\)T [ J{w,\);'{VMyt){^t,^t,dwt) 
Jo 

+ Va{y,){ildh) + Va{yt){ildht) + V%{yt){m^)dt]. (5.5) 
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We can also obtain higher order directional derivatives in a similar way, but we omit them. 

In a proof for the main theorem, we need to consider y^{w) = I{ew + 7, e^^^X), where 
7 G H is a fixed element and e G (0, 1]. This process satisfies the following Young SDE; 

dft = a{yl)edwt + (T{yl)dlt + h{yl)e^'^dt with yl = a. (5.6) 

When 7 = 0, we write for y'^. In that case, self-similarity of (wt) implies that the two 
processes (?/£i/Hj)o<t<i and {yf)o<t<i have the same law. 

In the next proposition we give estimates for the derivatives D'^y^. As we stated 
above, it is known that yx (and hence y^) is Dqo- In that sense, this proposition is not 
new. But, the estimate in powers of e in (15.71) may be new. Also, the proof is slightly 
different from the preceding papers, because Proposition 14.51 is used. 

Proposition 5.1 Take any 7 G H and fix it. Then, for any q G (l,c>o) and k = 
0, 1, 2, . . . , there exists a positive constant Cq^k such that 

¥.[\\D^y^\'' < Cq^ue^ for all e ^ (0,1] andT e [0,1]. (5.7) 



Proof. In this proof, an unimportant positive constant C may change from line to line. 
First, consider the case A; = 0. Since uj{s,t) = i\\w\\^_i^i^ + + l)(i — s) satisfies 

\{ewt + It) - {ews + ls)\ + \e'/''t - e'/^'sl < uj{s, tf'^, (0 < s < t < 1), 
we can use Proposition 13. H (ii) to obtain that 

Wt\ < WL^kid < \a\ + C{1 + \\wZ_^^, + ||7|IL, J (5.8) 

for some constant C = Ck- Fernique's theorem immediately implies (15. 7p for k = Q. 

Next let us consider the case k = 1. By slightly modifying (I5.2p - (]5.3p . we can easily 
see that ff'^ := Dtyf satisfies the following (l5.9l) -( l5rT0l) : 

dif-Vcx{yt){et'\ d{ewt+it))-Vh{yt){lf)e'/''dt = a{yt)edht with = 0, (5.9) 

or equivalently, 

= Jt / J,-'aiyt)edht, (5.10) 
Jo 

where J = J{ew + 7, e^^^ \). From this, we can easily see that 

a— hid < Cell Jllooll^. cr{y':)\\ Q— /lid II '''\\a~hld 

< Ce\\J\\oo\\J'^\\a-hid{l + W\\a~hid)\\h\\n (5.11) 

and, hence, 

II^^Tlk* <C£||J||oo||J~l a— hid 1 + II y'^Wa-hld)- 
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By Proposition 13.41 and Remark \3M L'^-norm of || J^^||a-/i/d is finite and bounded in e for 
any fixed g G (1, oo). (Just replace w and t in Proposition 13.41 and Remark [3l6] by ew + 'y 
and e^/^t, respectively.) Hence, using Holder's inequality, we obtain (15. 7p for A; = 1. 

We prove the case k = 2. Set DkDhy^ = 0^^'^ for simplicity. Then, in the same way 
as in (I5.4p - (l5.5p . we have 

If^^ = Jt f J;'{^Mytm\ et'\ d{ew, + 7*)) 
Jo 

+ Va{ytm\edh) + Va{yt){lf\sdh,) + V%ymfrQ'^)dt] . (5.12) 
From this, we have 

\\C'^'^\\(x-hid < c*!! j||oo||^"lQ-/i/d{||W(y^)|| 

a— hid a— hid \\C,T Wa-hldW^W + 'y\\a-hld 

+ ||Vo-(y;^)|| a— hid i\\C*\\a-hld\\£h\\a-hld + \\C'''\\a-hld\\£k\\a-hld) 

+ ||VV(y^)|U_,i,||f^i,_;,;rf||f^^|U_,;4 

< Ce'\\J\\UJ-'\\Lhidil + Ml^Hid + ht-hidnh\\n\\k\\n- (5.13) 
Here, we used (15.81) and (15.111) From Proposition 14.51 we see 

Wif'^Wn^^n* < Ce'\\J\\l\\J-ra^Ul + \Ml-hid + ht-hid)'- (5-14) 

Using the moment estimate for || again, we show (15.71) for k = 2. The higher 

order cases {k > 3) can be dealt with in a similar way. I 



Remark 5.2 We already have (i) Frechet smoothness of in the deterministic sense 
and (ii) L'^-estimates for derivatives as in this proposition. From these, we can easily 
verify that G Dqo as follows. (For simplicity of notations, we only consider the case 
7 = 0,e = l.y) By using Taylor expansion, we have 

yriw + rh) - yr{w) _ ^^^^^^^ ^ ^ ^^^^ _ ^^^2^^^^ ^ 



' JO 

for all w G W, h G Ti^ , and r G R. Note that the derivative D on the both sides of 
the above equation is in the deterministic sense. By Proposition \5. 1\ and Cameron- Martin 
formula, the right hand side is 0{r) as r ^ in L'^-norm for any q G (1, oo). This implies 
that yT G Dg i for any g G (1, oo) and the derivative Dyx in (deterministic) Frechet sense 
is also the derivative in the sense of Malliavin calculus. (See Proposition 4-21, /7P| / for 
instance.) The higher order derivatives can be dealt with in the same way. 

Now we show that, under the ellipticity condition (Al) for a (i.e., for Vi, . . . , V^), the 
Malliavin covariance matrix for 

^ (5.15) 
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is uniformly non-degenerate in the sense of Malliavin as e \ 0. Here, we set a' = for 
the solution of the following ODE; dcf)^ = a{(j)^)d'yt with 0q = a. 

(Al): The set of vectors {Vi(a), . . . , Vd{a)} linearly spans R". 

Nualart and Saussereau |T7j showed non-degenarcy of Malliavin covariance matrix 
of yT under (Al). Baudoin and Hairer [1], proved non- degeneracy under Hormander's 
hypoellipticity condition for vector fields {Vi, . . . , V^; Vq}. 

In the next proposition, we will prove uniform non- degeneracy of flS.lSp under (Al) 
by slightly modifying Baudoin-Hairer's argument. (The special case 7 = has already 
appeared in Baudoin and Ouyang [3J.) 

Proposition 5.3 Let = [y^'^, . . . ,y^'"') be the solution of ( 15. g|) and assume (Al). 
Then, {y{ — a')/e is uniformly non- degenerate in the sense of Malliavin as e \ 0. 



Proof. Let y = {yt) be the solution of (15. ip . In p. 388-389, [T], an explicit form of 
the Malliavin covariance matrix for yi is given. By replacing the vector fields [a; h] = 
[Vi, . . . , Vrf; Vo] with [ea; e^^^b] = [eVi, . . . , eVd, e^^^Vo], we can easily see that 



1 /•! 



Jo 



J{ew, e^'''X)Z^a{yl)a{yiyj{ew, s^/'' X)~^^*\u - vl^^'-^dudvJiew, e^/^A)^. 



Here, \t = t and A* denotes the transposed matrix of A. By shifting w ^-^ w + {•y /e), we 
have 

l((2?r,^^r^)«) = H{2H-l)J,C.r„ (5.16) 

where Jt = J{ew + 7, e^^^X)t as before and we set 

C= [' f J;^cy{yl)a{yiyj;^^*\s-t\'''-Hsdt. 
Jo Jo 

Since supo<£<i H-^i'^Hq < 00 for any q E (1, 00), it is sufficient to prove 

sup II I det C|^"'^||y < 00 for any 1 < g < 00. (5-17) 

0<£<1 

We will follow the argument in pp. 387-340, yy. In order to show (I5.17P above, it is 
sufficient to prove that, for any 1 < g < 00, there exists po{q), which is independent of e 
and satisfies that, 

sup /i((a, Ca) < p) < p'' for any p G (0,po(?)) and e G (0, 1]. (5.18) 

aGR",||a||=l 
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(For a proof, see Lemma 2.3.1, Nualart [IS]). As in [T], 



(a,Ca) = ^ r [\B.J:'VM)){^.Ji'VM))\s-tr-'dsdt 
Jo Jo 



d 

|2 



= 5^||(a,J.-V,(r))||J,. (5.19) 
i=i 

By a Norris-type lemma (Corollary 4.5, |1]), there exists < /3 < 1/2 such that for any 
7 < if — (1/2) and < p < 1, the following inequalities hold; 



/z((a,Ca) < p) < min u(||(a, J.-V,(r))||^ < P^^') 

p(||(a,Jr^V,<r))L. <P^/^)+/i(ll(a,JrV,(y^))||^_,^^>p-^/^)^ 
p(|(a,\/,(a))| < p^/^) +p(||(a, J.-V,(r))L_,,, > p-''% (5.20) 



< min 

l<i<d 



< min 

l<i<d 

Here, in the last inequality, we evaluated at t = and used 7 < a. Note that the set in 
the first term on the right hand side is already independent of e and non-random (i.e., 
either or the whole set W). 

Recall that, for any g, + lly^lla-hiJ — '^1 some constant Ci = Ci(g) 

which is independent of e. Then, using Chebyshev's inequality, we have 

p(ii(a,jrv,(r))|L-,,,>p~^/^) <C2P'' 

for some constant C2 = C2{q) which is independent of e. 

Let us consider the first term on the right hand side of fl5.20p . By (Al), there exists 
c' > such that a{a)a{ay > c'Id„ in the form sense. We have 

n max \{a,Vj{a))\^ > V \{si, Vj{a))\^ = {si,a{a)(x{a)* Si) >c' >0. 

l<j<d 



Hence, if p^/^ < \/dJn, then mini<j<dp(|(a, Vj(a)) | < p^/^) = and 

p((a,Ca)<p) <C2(g)p«. 

From this, we can easily see f lS.lSp holds with po{q) = C2{q + l)^^ A [d /nY^^ . This 
completes the proof. I 



Remark 5.4 Precisely speaking, another Hilhert space 'H is used in fl^ instead ofH. (See 
Appendix below for the definition ofH.) This causes no confusion, however, since the two 
Hilhert spaces are unitarily isometric. 
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6 On-diagonal short time asymptotics 



The aim of this section is to prove Theorem l2.1l namely, on-diagonal short time asymptotic 
expansion of the density of the solution of the Young SDE (12.11) (or equivalently (15. ip ) 
under the ellipticity assumption (Al). 

Let us consider the solution (yt) = {yt{,o)) of Young differential equation (12. ip with an 
initial condition yo = a E R" driven by fBm {wt) with H > 1/2. It is shown in [iTl [1] 
that, under (Al), the law of the solution has a smooth density p{t, a, a'), i.e.. 



P(yt(a) e A) = p{t, a, a')da' ( for any Borel set A cW 
J A 



)■ 



For t > 0, yt = yt{ci) is Dqo and non-degenerate in the sense of Malliavin. By the same 
argument as in Ikeda and Watanabe [8], we have the following expression; p{t,a,a') = 
E[6a'{yt{a))] = D_^(5a'(z/t(a)), 1)doo- By the self-similarity of fBm, (?/ei/fff)f>o and {yt)t>o 
have the same law, where is given by (15.60 with 7 = 0. From this, we see that 

p{e'/^,a,a')=E[6Ayl{a))]. 

The most important part of the proof is asymptotic expansion of in £ G (0, 1] in 
Doo-topology. For that purpose, we introduce the following index set for exponent of e. 
Set 

Al = {r^i + ^ I ni,n2 G N}. 

We denote by = < /ti < K2 < ■ ■ ■ the elements of Ai in increasing order. Several 
smallest elements are explicitly given as follows; 

1 1 2 

Kl = 1, ^2 = -^, ft^S = 2, K4 = l + — , K5 = 3A— 

Proposition 6.1 The family of Wiener functional yl {0 < e < 1) admits the following 
asymptotic expansion as e \0; 

~ a + e/i + e^^f^, + e^^f,, + ■■■ m Doo(R") 

for certain /^i, /^a, • • • G Doo(R"). 



Proof. For j = (ji, . . . , j„) G {0, 1, . . . , rf}™, we set |j| = m and 



H 



i{l<k<m\jk^ 0}. 



We denote by the totality of such j's with |j| = m and set X = W^^{Im- 

We will use the following convention. We set t = w^. Then, the ODE for (that is, 
(15. 6p with 7 = 0) reads; 

d 

dyt = e""VM)dw'l + ^ eVM)M with = a. (6.1) 
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It is easy to see that 

Jo Jo 





Here, ViVj denotes a vector field Vi (as a first order differential operator) acting on a 
R"- valued function Vj. 

Repeating the same argument for the last term on the right hand side of (16. 2p . we 
have 




Let us observe the first term. From basic properties of Young integral, we easily see 
that, for any m, the real-valued functional J^^ ■ ■ ■ /q™"^ '^''^tZ ' ' ' dw{ldw{l is in mth 
(inhomogeneous) Wiener chaos and hence it is in any Gq,k (1 < Q' < oo, k G N). 
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Next we consider the last term in (16.41) . We set 

rti rt„-i 

Qe{w)= / ■•■/ V,„---V,,V,MJdwl:---dwrjwil 
Jo Jo Jo 

and will prove = 0(1) as e \ in Dq^fc(R") for any 1 < g < oo, /c G N. (For simplicity 
of notation, we denote G = Vj^ ■ ■ ■ Vj^Vj^ and assume jj ^ for all i. The other case is 
actually easier.) 

Since is 0(1) in any L"? as e \ 0, Qe{w) is 0(1) in any L^, too. Now we 

estimate the derivatives. For h ^H, we have 

DnQsiw) = ■ ■■ VGiyl){D,yl)dwl: ■ ■ ■ dwUdwH 
Jo Jo 

+ E / / G{yl)dwl:---dhl[...dwll. 

Jo Jo 

Holder norms of y'^ and D^y^ were estimated in (I5.8p - (l5.1ip . From these, we see that 
WDQeWw = 0(1) in any L". 
Similarly, h,k E Ti, we have 

/■l rtn-i 

DkDnQeiw) = ■ ■■ VG{yl){DkD,yl)dw{: ■ --dwildwil 

Jo Jo 

rl rtn-i 

+ ■■■ / VGiyl){D,yl,D,yl)dwl:-- -dwildwil 

Jo Jo 

+ E / ■ ■ • / '^GiyD^Dr.yDdwl: ■ ■ ■ dk^ ■ ■ ■ rf< 

Jo Jo 

+ E / ■ ■ • / "^GiytJiD.yDdwi: ■ ■ ■ dh^ ■ ■ ■ rf< 

,_, Jo Jo 



1=1 

n „i Ml 



+ E / ■ ■ ■ / G{y!Jdwl: ■ ■ ■ dhl\ ■ ■ ■ dKZ ■ ■ ■ dwil. 

Holder norm of DkDhy'^ was estimated in (15.131) . Combined with Proposition 14. 5[ the 
above implies that ll-D^Qell-H*®-^* = 0(1) in any L^. Higher order derivatives can be done 
in the same way. 

Now we prove the proposition. In order to get the asymptotic expansion up to order 
Km (i.e., the remainder is of order k^+i), it is sufficient (i) to consider the expansion 
(16. 4p with n — 1 being the smallest integer which is not less than and (ii) to set 

/ ■■■ / dwi: ■ ■ ■ dw^dwii 

for all 1 < / < m. I 
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Before we prove on-diagonal short time kernel asymptotics, we define two more index 
sets for exponent of e. Set A2 = {k — 1 | k G Ai \ {0}}. Smallest elements of A2 are 

1 1,2. 
0, Jj-^. 1, -Jj. (3A-)-l,... 

Next we set A3 = {ai + 02 + ■ ■ ■ + | G N+ and ai, . . . , G A2}. In the sequel, 
{0 = z/q < z^i < z^2 < ■ ■ ■ } stands for all the elements of A3 in increasing order. 



Proof of Theorem \2.1\ First, note that 

Pie'/^a,a) = E[5,(yf (a))] = E[<5o (^^i^^)] = .-"E[5o(^i^^)]. 



By Proposition 15.31 {yl{a) — a)/e is uniformly non-degenerate. It admits asymptotic 
expansion in Doo(R'^) as in Proposition 16. II Then, by Theorem 14. 4[ the following asymp- 
totic expansion holds in D_oo as e \ 0; 

5o (Mi^in^) ^ 00 + 5-10^^ + e-^^^ + . . . as £ \ 0. 

By taking the generalized expectation and setting Cy^, = E[0y^], we have 

p{e^^^, a, a) ~ e"" (cq + c^.e'"^ + c^^^"^ H ) as e \ 0. 

Putting 6 = , we complete the proof of Theorem 12. 1[ I 



7 Taylor expansion of Ito map around a Cameron- 
Martin path 

In this section we prove an asymptotic expansion for y'^ = I(ew + 7,6^/^ A), which was 
defined in (15. 6p . The base point 7 G "H of the expansion is arbitrary, but fixed. First, we 
prove that admits the following expansion in C"~'^'^([0, 1]; R"); 

r ~0° + £"^0"^+e"2 0"2 + --- as£\0, (/€iGAi=N + ^N), (7.1) 

for some C°'^^^'^{[0, 1]; R")-valued Wiener functional (p^, 0*^^, 0"^^, . . .. Since the Ito map / 
in the sense of Young integral equation is smooth in Frechet sense (see [H]), this kind 
expansion holds in deterministic sense. In this paper, however, we need to prove this 
expansion in L'^-sense. 

Before we state the proposition precisely, we now give a heuristic argument to find an 
explicit form of 0^"". To find an ODE for 0*^ is easy. 

dyl = a{yl){edwt + d'yt) + b{yl)e^^^ dt with yo = a, 

d(j)^t = a{(f)'^t)d-ft with 0° = a. (7.2) 
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Set A0 := ~ (fp and put it in the above ODE for . Then we have 
d{(tP + A0) = a{(tP + /\<i)){edw + d-i) + h{(f + A(j))6^/^dt 

^ ^' ^ 1 ' ^ ^ 1 ' 

Assume A0 admits the asymptotic expansion (17.11) . Then, by putting it in the above 
equation and picking up the terms of order e**", we find an ODE for 0^^'". Note that 
0^^™ = for all m > 1. 

For Km = 1, 2, we can write down the ODEs explicitly as follows; 

d<Pl - Va(0°)(0^ d^t) = (y{<ft)dwu (7.3) 
dc^]/^ - Va(0°) (0^^, rf7*) = h{4)dt, (7.4) 

dcfl - Va{(tPMld^t) = Va{<ftWt,dwt) + ^VV(0?)(0^ c?7,). (7.5) 

Note that 0^/^ is independent of w, i.e, non-random with respect to /x. 
For Km > 2, 

^ VV(0° 



Va(0°)(0^,t/7*)=E E ^-^{<PT\---AT'=;dwt) 

k = l KijH hKi^,=Km-l 

+E E ^^(0r\...,0r^^7.) 



+E E ^^(0rs...,0r^cit). (7.6) 

fc=l Ki-^H h/tij.=Km-(l/-H') 

The summations in the first term on the right hand side is taken over all , . . . , Kj^. G 
Ai \ {0} such that + ■ ■ ■ + = — 1 hold. Kj^. = is not allowed. So, the sum 
is actually a finite sum. The second and the third terms should be understood in the 
same way. An important observation is that the right hand side of (17.61) does not involve 
0*^™, but only 0", 0^, . . . , 0'*'"-i. These ODEs have a rigorous meaning. So, we inductively 
define 0'^™ as a unique solution of (I7.3p - (l7.6p . 

If the right hand side of (l7.3p - (l7.6l) is denoted by dQI"", then 0"™ can be written 
explicitly as follows; 

^t- = Ai)t r J{i)T'dQt-, (7.7) 



where we set 7(7) = J(7, 0) = J{Qw + 7, O^^-'^A). See (I3.17P for the definition of J . 
Define the remainder term R'^rn+i,e 

We will estimate this remainder term in L'^-sense. 



26 



Proposition 7.1 For any m G N and q G (1, oo), ll^**"* Ho-ziid ^ L'^if^) (^iT'd 
E[\\R--^^'ro.-Hi,Y^' = 0{e''-^^) ase\0. 

Proof. From the expression (17. 7p and induction, it is easy to see that ||<^'*'"||a_/i/d G 
f^KqKooL'^ for any m. Let us consider i?^'*^ = A0 = y'^ — (jp = I{ew + 7, e^^^X) — 1(7, 0). 
Here, / stands for the Ito map and stands for one-dimensional constant path staying at 
0. 

Define u{s,t) = + ||7||a_wd + ~ ^) with a = 1/p. Since this control 

function satisfies 

\{ewt + 7t) - {ew, + 7.)l + l^'^^'t - e'/^'sl < 00(3, tf'^ 

\{{ewt + 7t) - {ew, + 7.)} - {it - 7.}| + le^'^'t - e^"" s\ < euj{s, tf'^ 

for all < s < t < 1 and e G [0, 1], we can use Proposition 13.21 to obtain that 

\rY - RY\ < eC{l + w(0, l))(P-^)/Pexp(Ca;(0, l))u{s, tf'' 

for some positive constant C. Since j9 < 2, we can use Fernique's theorem to obtain the 
desired estimate holds when k^+i = 1- 

Before we prove the higher order cases, let us observe the concrete expression for 
several i^'^'^+^'^'s. In the sequel, we write Km+i ='■ i^m+ for simplicity of notation. First 
we consider R^^'^ = R^/^'^ = y^ ~ (p^ — ecj?- . A straight forward computation yields; 

dR]^'' = e{a{yl) - a{<f,)}dwt 

+ [{a{yl) - a(0°)}rf7* - Va(0°)(e0i, ^7*)] + s^'^'biyDdt. (7.8) 

From this, we immediately have 

dR\^'' - Va{<p'l){Rl^'\d^t) = e{a{yt) - a{<j>'i)}dwt 

+ \f^ d9Vai<p', + 9RY){RY, R]'\ d^,) + e"%{yl)dt (=: dh]^^'). (7.9) 

Observe that, on the right hand side, there are only i?^'^, y^, 0°, 7, w, which are known 
quantities, but no Since R}j^'^ = J{i)t Jq J{l)7^dLl~^''^ as before, it suffices to 

show that \\L^+''\\a-hid = 0{e^/^) for any L'?. 

Since \\e^/^ J^biyDdtW^-hid < Ce^/^\\y'\\a-hid, the third term of L^+'' is 0{e^/^) in 
any L?. Similarly, e|| J^iaiyl) - a{(j)°)}dwt\\a-hid < C€\\R^'^\\a-hid\\u)\\a-hid, the first term 
ofLi+'^isO(£^)inanyL«. For any ^, ||Va(0O + ^/?i'^)|U_/,zrf < C(||0O|U_;,zrf+pi'^|U_/,zrf). 
Hence, we have 




1 



^0 



deVa{<p'i + dRl''){Rl'',Rl'',d-ft)\\^^Hid<C{\\<p''\\^.Md+\\R'''^ 
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hid- 



We see from the above inequality that the second term of L^"*"'^ is 0(£^) in any L'^ and 
hence Hi^^^'^Ha-ftid = 0{e^/^) in any L^. Thus, we have obtained the desired estimate for 

The estimate for R^^l^)^'^ = — cj)^ — ecj)^ —e^^^ (j)^/^ can easily be obtained as follows. 
We can immediately see from (17. 5p and (17.91) that 



2 - „ 

+ e'^'mi) - (=: rfLf /^)+'^). (7.10) 

Notice that we have essentially shown that = 0{e^) in any L^. Thus, we 

have obtained the desired estimate for = i^^.e^ 

Next, we will estimate = _ 0O _ ^^i _ ^i/h^i/h _ ^2^2 _ p^^^ (^^^ (^^^ ^^^^ 

(17. 8p . we see that 

dR',+'' = [{aivt) - ai<p',)}edwt - Vai<p',){e^ledwt)] 

- ^VV(0°)(e0,\e0i,rf7*) - (7.11) 

The second term on the right hand side is equal to 

Va(0?)(^ii?r^^7*) + lv'a{^',){R]'',Rl'',dj,) 

+ I' ^V V(0° + ^^i?^^)(i?,^-, i?,^^ i?,^^ ^7.). 

Hence, (17.111) is equivalent to the following; 

cii?^'^ - Va(0°)(di??+^rf7i) = [Wift) - (r{4>'t)}edwt - Va{^',){e(Pledwt)] 

+ ^[vv(</.?)(i^^^i^^^rf7t) - vv(0°) (50^501,^7*)] 

+ I' ^V V(0° + ^^i?^^)(i?,^^ i?,^^ i?,^-, rf7*) 

+ ^'/""{KyD - bi<l>'t)}dt (=: ciL^'^). (7.12) 

Then, 4+'^ = J(7)t/o^ J(7)r'rf^^'^ 

Let us observe the right hand side of (17.121) . There are no 
assumption of induction, we may only use the relation R'^'^ = = — (p^ — ecj)^ — 

£^/^(f)^/^ and the estimates of R'^'^ for k = 1, 1/H,2 (and of ^'"s). In the same way as 
above, by using the Taylor expansion, we can prove that ||i^^''''^||a-ftid = 0{£^^^^^^^) in 
any L'^. Cancellation of the terms of order < 2 on the right hand side is no mystery 
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because of the way are defined. Tlius, we have obtained the desired estimate for 

Higher order remainder terms can be dealt with in a similar way. We give a sketch of 
proof. There exists 

such that c/i^^^'"+''^-Va(0°)(rf/^^'''"+''^^^7^) = dLt"'+'''. Due to cancellation \\L''"^+^''\\^.hid = 
0{e'^"'+^) holds in any L'^. This proves the assertion. I 

The next proposition shows that, when evaluated at t = 1, Eq. (17. ip gives an asymp- 
totic expansion in Doo(R"). 

Proposition 7.2 We have the following asymptotic expansion in Doo(R'^). 

+ + ase\0. (7.13) 

Here, = < ni < K2 < ■ ■ ■ c^re all the elements o/ Ai = N + -^N in increasing order. 

Proof. By using induction and basic properties of Young integral, we can easily see that 
(pl"^ is in [Km]-th inhomogeneous Wiener chaos for each t and m. In particular, 0^™ G Dgo- 
If A; > + 1, then D''RI^+''^ = D^yl From Proposition EH this is 0{e^), and hence 
Q(^^Km+i^ in any L^. The Meyer equivalence implies that R'^"'+^''^ is 0{e'^"'+^) in D^^^ for 
any q and sufficiently large k. Since ^-norm is increasing in k, the proof is completed. 
I 

We now recall the following Taylor expansion of Ito map around 7 in the deterministic 
sense. 

Lemma 7.3 (i) For each m, there exists c = c(Km) such that 

H^-\\a-hid < c(l + M^^Mdr- for all w e Co"-'^'"([0, 1], R^). 
(ii) For each m and r > 0, there exists d = c'{K,m,f) such that 

\\R''"'+'''\\a-hid < c{e + \\ew\\a-hidT"'+\ if \\ew\\a-hid < r. 

Proof. This is immediate since = I{ew + 'y,e^^^\) and Ito map I is Frechet smooth 
by Li and Lyons's result [H] . It is also possible to prove this lemma by using the explicit 
expression of and mathematical induction as in the proof of Proposition l7. II above. 

I 
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8 OfF-diagonal short time asymptotics 



In this section we prove the short time asymptotics of kernel function pt{a, a') when a 7^ a'. 
We basically follow Watanabe [20]. In this paper, however, we can localize around the 
energy minimizing path in the abstract Wiener space since Ito map is continuous in our 
setting. This makes the proof slightly simpler. 

8.1 Localization around energy minimizing path 

For 7 G 7{, let 0° = be a unique solution of (17.21) . which starts at a G R". Set, for 

a 7^ a', 

K' = {7 G ^ I 0?(7) = a'}. 

We only consider the case that K^' is not empty. For example, if (Al) is satisfied for any 
a, then K^' is not empty for any a'. From the Schilder-type large deviation theory, it is 
easy to see that inf{||7||^ | 7 G K^'} = min{||7||-^ | 7 G K^'}. 

We continue to assume (Al). Now we introduce another assumption; 

(A2): 7 G K^' which minimizes "H-norm exists uniquely. 

In the sequel, 7 denotes the minimizer in Assumption (A2) and we use the results of the 
previous section for this 7. 

Note that (i) the mapping 7 G W 1— )■ 0?(7) G R*^ is Frechet differentiable and 

(ii) its Jacobian is a surjective linear mapping from "H to R'' for any 7, because there 
exists a positive constant c = 0(7) such that 

W'^(7),I^0?'^(7))kO >c-Id„. (8.1) 

\ / l<i,j<n 

This can be shown in the same way as in the proof of Proposition 15. 3[ (Actually, it is 
easier since 7 is non-random and fixed here.) 

Therefore, by the Lagrange multiplier method, there exists u = {ui, . . . ,i?n) £ R" 
uniquely such that the map 

nxR"3 (7, z.) ^ ^W^Wl - {u, 0?(7) - a')R. G R (8.2) 

attains extremum at (7, u). By differentiating in the direction of G "H, we have 

(7, k)n = {i^, I^fc0?(7))R" = ^(7)1 m^M<Ptmdkt)^^. (8.3) 

Jo 

Here, the definition of ^(7) was given just below (17. 7p and the integral on the right hand 
side is Young integral. Hence, (7, ■ extends to a continuous linear functional on W. 
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Let us introduce Besov-type norms. In the context of Malhavin calculus, these norms 
are often more useful than Holder norms and p- variation norms since (a power of) these 
norms become Doo-functionals. For m > 0, < < 1, and x G Co([0, 1], R'^), we set 



\X\\m,e-B 

l0<s<t<l 

Ym.e-B, 



dsdt 



and Co™'''"^([0, 1], R"') = {x G Co([0, 1], R'^) | \x\mfi-B < oo}. It is known that \\x\\e^hid < 
c\\x\\m,e-B for some constant c = Cm,e > 0. Hence, C™'^~^([0, 1],R<^) C ^^^"^'^([0, 1],R^). 

Let (wt) be fBm with Hurst parameter H G (1/2, 1) and let a(= 1/p) < H as before. 
Since K[\wt — w^p] = d\t — s]"^^ , we can easily see < oo if m > 1/{H — a). 

Therefore, the law of fBm, /i = /x^, is supported in C™'"-^([0, 1], R'^) if m > 1/{H — a). 
Set Wb to be the closure of Cameron-Martin space H = in C'°"''([0, 1], R'^). Then, 
{yVB,'H,fi) is also an abstract Wiener space. 

Now we recall Schilder-type large deviation principle for scaled Gaussian measures. 
For e > 0, let fis be the law of the law of the process {ewt)o<t<i- This is a measure 
on Wb- Set X(w) = ||w||^/2 (if w ETi) and I{w) = oo (otherwise). It is well-known 
that X : Wb — ?■ [0, oo] is lower semicontinuous and that X is good, i.e., the level set 
{w I X(w) < r} is compact in Wb for any r G [0, oo). 

The family {/ie}e>o satisfies large deviation principle as £ \ with a good rate 
function X, that is, for any measurable set A C Wb 

— inf X(w) < liminf log/ie(A°) < lim sup£2log;Ue(A) < - inf I{w). (8.4) 

w€A° e\0 e\^0 w6A 

Next, set fie = l^e ® K^/H\, where A is a one- dimensional path defined by = t 
and ® stands for the product of probability measures. In other words, fie is the law of 
the {d + l)-dimensional process {ewt-, o<t<i under /x. This measure is supported on 
Wb © R(A) C C™'""^([0, 1]; R'^+i). Define X(m;; /) = ||w||^/2 (if w G H and k = 0) and 
X(w, /) = oo (otherwise). Here, / is a one- dimensional path. 

From (18. 4p we can easily show that {fte}e>o satisfies large deviation principle as e \ 
with a good rate function X, that is, for any measurable set A C Wb © R(A), 

— inf i{w) < liminf log/ie(y4°) < lim sup log /ig (A) < — mf_X{w). (8.5) 

w€A° e\0 £\^0 w£A 

We will use (18. 5p in Lemma \87i\ below to show that we may concentrate on a neighborhood 
of the minimizer 7 contributes to the asymptotic expansion. 

From now on, we will fix an even integer m > such that m > 1/{H — a). Then, it 
is easy to check ^ Dqo. In fact, this functional is an element of mth inhomoge- 

neous Wiener chaos, i.e., D'^~^^\\w\\'^ = 0. 

Now we introduce a cut-off function. Let : R — [0, 1] be a smooth function such 
that ipiu) = 1 if \u\ < 1/2 and ipi^u) = if \u\ > 1. For each r] > and e > 0, we set 
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We can easily see that Xvi^y ' ) ^ ^oo- Shifting by •j/e, we have 

It is easy to see from Taylor expansion for that, for any 77 > and any M G N, the 
following asymptotics holds; 

X^ie, w + ^) = l + O(£*0 in Doo as e \ 0. (8.6) 

The following lemma states that only the paths sufficiently close to the minimizer 7 
contribute to the asymptotics. 

Lemma 8.1 Assume (Al) and (A2). Then, for any r] > 0, there exists c = c,, > such 
that 

II 1 1 2 

< E[(l - w)) ■ Sa'ivl)] = 0(exp{-^^^^}) ase\0. 



Proof. We take rj' > arbitrarily and we will fix it for a while. It is obvious that 



0<E[{l-Xr,{e,w))-6a'{yl)] =E {1 - Xr,{e , w))iP 



m - a 

rj''- 



l\2 



5a' {yl) 



1.7) 



Set g{u) = u V for u G R. Then, in the sense of distributional derivative, g"{u) = 5o. 
Take a bounded continuous function C : — )■ R such that C(mi,...,m„) = g{ui — 
a[)g{u2 — a'g) ■ ■ ■ g{un — a'^) if \u — a'\ < 2rj' . Then, the right hand side of f l8.7p is equal to 



E 



(l-^)( — 1|£^-7""^ 



m,a—B 



\yl - a ' 



l\2 



{dl---dlC){yZ 



Now, we use integration by parts on the abstract Wiener space as in f l4.ip - fl4.2p to see 
that (18. Sp is equal to a finite sum of the following form; 



5^E[F,-fc(e,i/;)-(l-^)(^)(-l^||e^-7| 



m 
m,a—B 



(k)f \yi-(^' 



i\2 



12 



.9) 



Here, Fj^is{e,w) is a polynomial in components of (i) and its derivatives, (ii) \\ew — 
iWma-B ^nd its derivatives, (iii) r(e), which is Malliavin covariance matrix of yf, and its 
derivatives, and (iv) 7(5) = T{e)~^. Note that the derivatives of 7(5) do not appear. 

From Proposition 15.31 there exists r' > such that |7*-'(£)| = 0{e~^ ) in L'' as £ \ 
for all 1 < g < 00. (Recall a well-known formula to obtain the inverse matrix with 
the adjugate matrix of A divided by detA.) Therefore, there exists r > such that 
|Fj,fe(£)| = 0{e''') in L"? as e \ for all 1 < g < 00. 
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By Holder's inequality, (18. 9p is dominated by 



c 



5^e[|(i-^)(^-)(- 



\ew — 7| 



m,a—B j 



/|2 



^.10) 



Here, l/g + l/g' = l and c = c{q,q',r],T]') is a positive constant, which may change from 
line to line. 

Since we may let q' \l after taking limsup, we obtain the following; 



limsup£^logE[(l - Xvi^,'^)) ' ^a'ivl 

e\0 



< limsupe log/i 

e\0 



> 



1] 



Iz/i - « I < ^7 



lim sup log /i^ {(ti?, /) G © R(A) 



< - mi 



(8.11) 



Here, / denotes the Ito map corresponding to ODE (15.11) and we have used the large 
deviation for the last inequality. (Note that continuity of Ito map is used.) Recall that 
00 [^] = 0) is given by ODE (17:21) . 

Now let rj' tend to 0. As rj' decreases, the right hand side of (18. lip decreases. The 
proof is finished if the limit is strictly smaller than — 1|7||^/2. Assume otherwise. Then, 
there exists {7fe}^i C % such that 



|7.-7C,„-B> 



1 



^°[7fc]i - a'l < T, and, liminf( 



> 



mi 



In particular, {7^} is bounded in T-L and, hence, precompact in W^. Let 700 be any limit 
point. For simplicity, a subsequence that converges to 700 is again denoted by {7^}. Since 
7 1-^ 0°[7]i is continuous with respect to the topology of Wb, we see that 0°[7oo]i = cl' 
holds. Also, we have ||7oo —'iWma-B — v"^/'^- So, 700 7^ 7. From the lower semicontinuity 
of the rate function, we see that 700 € Ti and ||7oo|||{/2 < ||7|||^/2. This clearly contradicts 
Assumption (A2). I 



8.2 Integrability lemmas 

In this subsection, we prove a few lemmas for integrability of Wiener functionals of expo- 
nential type which will be used in the short time asymptotic expansion. 
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Throughout this subsection we assume (A2). Let 7 be as in (A2) and let 0"^ and 
^Kj+,e _ = 0, 1, 2, . . .) be as in Section [7] with 7 = 7. First we consider 



= e'^^- V" + £"^"^0"' + ■ ■ • . 
Here, K4 = 1 + iMB) and = 3 A (2//7). 

Lemma 8.2 Assume (A2). For any M > 0, there exists rj > such that 
sup E[exp(M(z7,/??+'^)/£2)j^ll^^ll <oo. 

0<e<l 

Proof. By Lemma FTSl if ||eu'||Q,_/ii(i < 1, then there exists a constant Ci, C2 > such that 
Hence, if < < 1, then 

Recall that, by Fernique's theorem, there exists a positive constant /3 > such that 
E[exp(/3(1 + \\w\\m,a-B)'^)] < oo. Take < 77 < 1 so that M\u\c2{2rjY^/"^-'^ < (3. Then, 
we see that 

sup E[exp{M{i^,Rl+'')/e^)l[\\eyj\u^_s<r,}] < 00. 

0<e<ri 

Note that, if ||£w||m,a-B < ''7 and 77 < £ < 1, then H-R^+'^Ha^/jirf/e^ is bounded. This 
completes the proof. I 

Next we consider 

^ = - 00 - £01) = + ^102 ^ . . . _ 

Lemma 8.3 Assume (A2). For any M > 0, there exists rj > such that 
sup E[exp(M||i?i+'^||t,w/^')/{||.»IU,._,<,}] <oo. 

0<e<l 



Proof. By Lemma [7.31 if H^wHa^/iifi < 1, then there exists a constant ci > such that 

ll^^^'lU-hid < ci{e + ||ew|U-ftZd)^^'^ < C2(£: + \\ew\\m,a-By^^- 
Hence, if ||£w||m,a-B < < 1, then 

\\R'^''\\l-,M/e^ < C2(l + \HU^.nne + vY^/"^-\ 
Then, we can prove the lemma in the same way as in Lemma 18. 2[ I 
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From now on we assume (Al) and (A2). In addition, we introduce the following 
assumption; 

(A3)': E[exp((z/, | 0} = 0] < cx). 

For all I < J < n, 0}'-' E C n*. When we regard ({){■' as an element of H 
by Riesz isometry, we write "0]^'"' El-id W^. We have an orthogonal decomposition 
T-i = kei © (ker^j^)-*-. We denote by vr the orthogonal projection from T-L onto ker0}. 
Note that (ker^J)-*- is an ra- dimensional linear subspace spanned by {V}'^; • • • ? "^i'"}- 
Since dim(ker ^j;)-*- < oo, the abstract Wiener space splits into two; Wb = ker^];" ^© 

(ker0j)^. The projection vr naturally extends to the one from Wb onto ker^j" ^, 
which is again denoted by the same symbol. There exist Gaussian measures ni and fi2 such 

that (ker0}" ^,ker0j,/ii) and {(kei (pl)-^ , (kei , H2) are abstract Wiener spaces. 
Naturally, fii = ix^ii, fi2 = ''^^^J' and /i = /ii x (the product measure). One may think 
Hi is the definition of the conditional measure P[- 1 0j = 0] in (A3)' above. 
Therefore, (A3)' is equivalent to the following; 

E[exp((z/,0? ott))] < 00. (8.12) 

Set 

+ ^^77)i^'^(7)r'VM0?)(0K^),0K^'),^7*), (8.13) 

where (pxiw) = J{^)t Jq J{'l)t'^(^{4>T)dwt- Then, ■j/' is a bounded bilinear mapping on Ws 
and so is 'iI}{71-,tt-). Clearly, ip{w,w) = </'f(w) and ifjinw^nw) = (pKnw). By Goodman's 
theorem (see Theorem 4.6, p. 83, P), restricted onH x Ti, (z/, ■iP{tt-, it-)) is of trace class 
and, in particular, Hilbert-Schmidt. The corresponding trace class operator on "H and 
corresponding element of the second Wiener chaos are denoted by A and S^, respectively. 
Then, {u , (j)l{7iw)) = S^(w)+Tr(A). Hence, (18.121) is equivalent toE[exp(S^)] < 00, which 
in turn is equivalent to sup Spec(A) < 1/2. Since the inequality is strict, there exists r > 1 
such that supSpec(ry4) < 1/2. This implies E[exp(Sj.yi)] = E[exp(rSA)] < 00. Summing 
it up, we have seen that (A3)' is equivalent to the following; 

E[exp(r(z/, (pf o tt))] < 00 for some r > 1. (8-14) 

Let us check here that (A3) and (A3)' are equivalent under (Al), (A2). 

Proposition 8.4 Under (Al) and (A2), the two conditions (A3) and (A3)' are equiv- 
alent. 
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Proof. As is explaind above, (A3)' is equivalent to supSpec(A) < 1/2. Keep in mind 
that the only accumulation point of Spec(y4) is 0, since A is of trace class. Let (— ^o) 3 
u I— 7- f{u) G K^' be a smooth curve in K^' such that /(O) = 7 and /'(O) 7^ as in (A3). 
Then, a straight forward calculation shows that 



'2 



\\f{u)\\l (\\f{u)\\ 



\H_ _ 

u=o 2 du^ 



u=0 



(z^,0;(/„)-a') 



= Wf'ml, - (z^,DV?(7)(vr/' (0),7r/'(0))> = ||r(0)||^ - (z^, ^(7r/'(0), 7r/'(0))>, (8.15) 

where we used f l8.2p -( l873l) and the fact that /'(O) is tangent to the submanifold K^' . Now, 
it is easy to see that supSpec(A) < 1/2 is equivalent to that right hand side of (18.151) is 
strictly larger than ||/'(0)||y2. I 

The following is a key technical lemma. It states that, restricted on a sufficiently small 
subset, exp{{h', R{'^) /e'^) G Ui<g<oo-Z^^ uniformly in e. 

Lemma 8.5 Assume (Al), (A2) and (A3). Then, there exists ri > 1 and rj > such 
that 

sup E[exp(ri(z/,i?i'")/e2)/|||,^ll^ <^|/^,^i../^l J < 00 

0<e<l II 1 / I- / J 

for any 771 > 0. 



Proof. By Lemma 18.21 and the relation Rl'^/e"^ = (pf + Ri'^'^/e'^, it is sufficient to show 
that 

sup E[exp(ri(z/,02))/{|,,^ll^ _^<^}/ 1,. ] < 00. (8.16) 

0<e<l u 1 / I- / J 

We give an explicit expression for the projection vr. Set Cjj' = 0}''' )-h* and 

^ ~ i^jj')i<j,j'<n ^ GL(?T,,R). The components of its inverse is denoted by C^^ = 
{Djji)i<jji<n- By straight forward calculation, tt : Ti ^ ker0j is given by 

nh = h-J2n4<l)l'',h)nD,,, 

hi' 

From this, it is easy see that vr : Wb — ?■ ker (f)\ is given by 

Tiw = w -^(l)]:\w)Djj, ■ ^Y- (8-17) 

Then, we have 

(j)l[w) = iIj{w, w) = (pliiTw) + 2 ^ (j)l'\w)Djj, ■ ij{w, 

jd' 

+ J2 <Pi'{w)(pi\w)D,,,Dkk'-i^{^<Pl''\^<l>Y')=-.Ji + J2 + J3. (8.18) 

,k,k' 
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Exponential integrability of the first term Ji on the right hand side of (18.181) is given 
in fl8.14p . So, we estimate the second term J2. Since e(j)\{w) = R\^''^{w) — RI'^{w), 



\<l)l'\w)^{wM'')\<ci 
for any c' > 0. 



e 



l,e/ 



\W\\m,a-B 



\W\ 



m,a—B 



Ac 



.12 





RY{w) 


} +Ci 


e 



^ \\'m,a—B 



Set C2 = 2cin^ supj J/ |-Djj'| and let M > 0. Then, by Holder's inequality, 



X 



/(4c')l 1/3^ r 3Afc2r,i|| 

II m,a — S 1 

1/3 



For any M > and r/i > 0, the third factor is integrable. If c' is chosen sufficiently large, 
then the second factor is also integrable by Fernique's theorem. By Lemma 18. 2^ there 
exists ?7 > such that sup^ of the first factor is finite and, hence. 



sup E[e^^l-^2l/{ 



<ri}h\RY/e\<m]\ < 



(8.19) 



0<e<l 

Since {w)4)\^ {w) = e^^{Rl'^''^ {wY — Ri'^ {wy}(f>l''^{w), we can deal with J3 in the 
same way. For any M > and ?7i > 0, there exists rj > such that 



sup Efe^^l-^^IJi 



0<e<l 



(8.20) 



Let r > 1 be as in (IHTTij) . Set n = (l+r)/2 > 1, g = 2r/(l+r) > 1, and 1/g+l/g' = 1. 
Then, from Holder's inequality and (18.141) . (I8.18p -( 18.20p . we can easily see that 



E [exp (n (Z/, (pl))l{\\ew\U,^.s<v}h\R'^'ye\<vi}\ 



< 



E[exp(r(z7,0?o7r))]'/'^J]E[e2'''^il^ll-^^l/. 



•j\\m,o.-B<v}^{\Rl'ye\<rii}\ 



1/(V) 



i=l 



From this, (I8.16P is immediate. This completes the proof. I 



8.3 Proof of off-diagonal short time asymptotics 

In this subsection we prove Theorem 12. 2[ namely, off-diagonal short time asymptotics of 
the density of the solution {yt) = {yt{a)) of Young ODE (15.11) driven by fBm (wt) with 
1/2 < if < 1 under Assumptions (A1)-(A3). 

First, let us calculate the kernel p{t, a, a'). Take r/ > as in Lemma [8.51 Then, we see 

pie'/'',a,a')=E[6aiyl)] 
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As we have shown in Lemma I8.H the second term I2 on the right hand side does not 
contribute to the asymptotic expansion. So, we have only to calculate the first term Ji. 
By Cameron-Martin formula, 



h = E[exp(-M| _ ^_(^^^^))sAyl)Xv{e,w + ^)]. 
Recall that (7, w) = {u, (p\{w)) for all w. Hence, noting that 0^/^ is non-random, we have 



h = exp(-M|)E[exp(-l(z7, 0}))5,.(a' + e<p\ + e^/^^;/^ + Rl^Xvie, 



lexpf-M + M^ 

X E[exp((z7, Rl^')/e')So{^\ + ^(V^)-^!/'' + e-'Rlnx,{e, w + ^)] 



e ' 



-^exp(-^ + ^)E[F(.,u,M„(^)]. 

where 

F{e,w)=exp{e-^u,R'/))Xr,ie,w + ^),p(\y^') (8.21) 

for any positive constant rji. It is easy to see that (i) Xvi^^w + 7/^) and its derivatives 
vanish outside {||£:u^||m,Q-B < v} (ii) V'(^r^|(^i ~ '^')/^| ) its derivatives vanish 
outside < r/i}. Hence, by Lemma [83| F{e,w) G Doo and F{e,w) = 0(1) with 

respect to that topology. Since So{{yl — a')/e) admits an asymptotic expansion in D_oo, 
the problem reduces to whether F{e,w) admits an asymptotic expansion in Dqo- 

Lemma 8.6 Assume (A1)-(A3). For any M G N, we have 

E[F(.,z.)5o(^^)] =E[F(.,u;)^(|</.l/r^i|^)5o(^^)] + 0(.^) 

as e \0. 



Proof. By using Taylor expansion for ip, we see that, for given M, there exist m G N and 
Gj{e,w) G Doo (1 < j < ^) such that 

H?l^r) - Hi^r) + i:^«(i^i>,(..u.) + o(e-) (8.22) 
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in Doo as e \ 0. Gj{6,w) = 0(1), but its explicit form is not important. Note that 
ip^^W(f)\/r]i\'^)T{(f)\) = if J > 1 and supp(r) C {a G R" | \a\ < 1/2}. 

By Theorem 14.41 and Proposition 15.31 So{{yl — o.')/^) admits an asymptotic expansion 
in D_oo as follows. As before, we set {0 = i^o < J^i < < • • • } to be all the elements of 
A3 in increasing order. For given M, let Z G N be the smallest integer such that M < z/;+i. 
Then, for some G D_oo (1 ^ J ^ 0) holds that 

- a')/e) = 5^{<P\) + e''^^,, + ■ ■ ■ + e^^^,^ + 0{e''^+-) (8.23) 

in D_oo as £ \ 0. Here, '^i,. is a finite linear combination of terms of the form 

9"5o(0i) X polynomial of the components of (pl" 's}. 

Hence, r]i\^)'^ uj vanish for all j, j'. 

Now, using f l8.22p and f l8.23p . we prove the lemma. 



= E[F(£,«;)^(l|^l^|')5o((yf - a')/e)\ 
= E[F(£,^)^(|0}^in5o((yf - a')/5)] 



+ E[F(e,^)(j]^(^)(Pf)G,(5,^))5o((yf-a')/^)] +O(e"0 
¥.[F{e,wm<P\h,\^)5,m-^')/^)\ 

rn ,1 

+ E[F(e,«;)(^V'^^)(|n|2)c',.(^,^))(5^(0i) + ... + ^^.$^J]+O(e 



Thus, we have shown the lemma. I 

Set A'2 = {/€ - 2 I K G Ai \ {0, 1, l/H}} = {0 < ff-^ - 1 < (3 A 2H-^) - 2 <■■■}. 
Next we set A3 = {ai + 02 + ■ ■ ■ + | m & N+ and ai, . . . , G Ag}. In the following 
lemma, {0 = po < Pi < P2 < ■ ■ ■ } stands for all the elements of A3 in increasing order. 

Lemma 8.7 Assume (A1)-(A3) and let F{e,w) G Dqo in Ii8.21\} . Then, for every 
/c = 1, 2, 3, . . 

Fie,wm<Pl{w)/r^,\') 

= exp((z7, 0?(«;))V^(|0i(u;)/r7in2{l + e'^^p^iw) + ■■■ + e'"'^,,iw)} + F^+iie, w), 
where Fk{e,w) G Dqo satisfies that 

Fk+i{e, w)T{(j)l) = 0{eP^+') m D_oo ase\0 
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for any T E 5'(R") with supp(T) C {a G | \a\ < 771/2}. Moreover, 7p^, G Dqo (j = 
1,2,.. .) are determined by the following formal expansion 1^^ = + 1 ); 

E = E ^{^"~'(^. + ^""'(^' + ■ ■ ■ } 

— ^ m! — ^ ml I ) 

m=0 m,=0 

Proof. Let ri > 1 be as in Lemma [8. 5 [ First we show that, for any ?7i > 0, 

E[exp(ri(i?,0?))/{l^i|<^^}] <oo. (8.24) 

We can choose a subsequence {sk] such that, as A; — )■ oo, \ and rI"^*" /sk — 0i a.s. 
To prove (18.241) . we apply Fatou's lemma to (I8.16P with rji replaced by 2rii. 

> E[exp(ri(z/,02)) liminf ^ E[exp(ri(z/, 02))/||^i|<^^}] . 

From (I8.24p . it is easy to check that exp[{iJ,(f)l{w)))ip{\(j)\{w)/r]i\'^) G Doo- 

Now we expand exp{{u, R{^) /e'^) = exp{{iy, (pl{w))) exp{{iy, R^^'^) /e'^) in e. Set Qi+i : 
R ^ R by 

We will prove that, for sufficiently large / G N, as e \ 0, 

e<^'<^?>gm((^,^i'''')/^')x.(^,«^ + ^)^(l0}H/^i|') = O(5'''=+^) in Doo. (8.25) 

Note that Xvi^^w + J) = 0(1) in Doo as £ \ by ([81]). By Proposition [TiSl Rl^'^/e^ = 
0(e(V^^)-i) in Doo. So, if / + 1 > pk+i/{{l/H) - 1}, then {{u, i??+'')/e2)'+i = 0{6p^+^) in 
Doo- Therefore, in order to verify (18.251) . it is sufficient to show that, as e \ 0, 

f\l _ gy^{p,<Pl+eRl+Ve')dQ . ^^(^e,w + ^)ij{\^l{w)/r],\^) = 0(1) in Doo- (8.26) 
Jo ^ 

To verify the integrability of this Wiener functional, note that e^" < 1 + e" for all m G R 
and < 6' < 1. This implies that the first factor on the left hand side of (18.261) is 
dominated by e^"'^^^ + e^'^'^i' ^Z*^^. From Lemma 1831 and (18.241) . we see that the left hand 
side of (I8.26P is 0(1) in any (1 < g < oo). In the same way, the Malliavin derivatives 
of the left hand side of (18.261) are 0(1) in any L^. 
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It is easy to see that, as e \ 0, 



E ^ ^^P'''^' = 1 + e'^l,, + ■■■ + e'^l,, + Oie^'^-^) in D... (8.27) 



fe=0 



From this and (18. 6p . we see that 

F{e,wm4>\{w)/r^^\^) 



= exp((j^,0?(w;))V^(|0i(t^)/r7i|>(^-2 
Using (18.221) . we finish the proof. I 



1 yl-a' ^ 



Proof of Theorem \2.2\ Here we prove our main theorem in this paper. We set 

A4 = As + A3 = {^^ + P k e As, p G A^;}. 

We denote by {0 = Aq < Ai < A2 < ■ ■ ■ } all the elements of A4 in increasing order. 
There is no mystery why this index set appears in the short time expansion of the kernel 
because, very formally speaking, the problem reduces to finding asymptotic behavior of 
E[exp{{u,Rl'')/e^) ■ 6o{Rl'ye)], as we have seen. Now, by (Km . LemmaEZl and (Km . 
we can easily prove Theorem 12.21 I 



9 Appendix 

In this appendix, we show that the result in Baudoin and Ouyang [2] is included in our 
main result. Let us recall the assumptions in p]. They set n = c? and assume (Al) for 
any starting point a G R*^ and, moreover, the following assumption (H): 

(H): There exist smooth and bounded real- valued functions ulj such that 

d 

colj = -uii and [Vi, Vj] = ^ u^jVi for all l<i,j,l< d. 

i=l 

Note that Vq does not appear in this condition. Under (Al) for any a, a{a)a{a)* is a 
d X d positive symmetric matrix, where cr(a) = [Vi(a), . . . , Vrf(a)] as before. As a result, 
a Riemannian metric tensor {gij{a))i<ij<d is defined on R'^ by g^^{a) = [a (a) cr (a) *]*•'. In 
terms of Reimannian geometry, (H) is equivalent to the condition that V^Y = [X, Y] 
for all smooth vector fields X, Y, where V"^*" is the Levi-Civita connection for this metric. 
From this, one can guess that this assumption is not very mild. 
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In [2], they proved short time kernel asymptotics under these assumptions when a 
and a' are sufficiently near. In the following proposition, we will prove that the main 
theorem (Thorem 1.2) in [2] is included in our main result as a special case. (In [2], unlike 
in this paper, the drift term Vq is zero. But, it does not make much difference in this 
proposition.) 

Proposition 9.1 Set n = d and assume (Al) for any a G R"' and (H). Then, if a and 
a' are sufficiently near, the conditions (A2) and (A3) are satisfied. 

The rest of the appendix is devoted for showing this proposition. So, we will keep 
assumption of Proposition 19.11 We basically follow Baudoin-Ouyang's argument, because 
it is in a sense shown in pp. 773-777, [2]. The differences are that (i) we use the Lagrange 
multiplier method (18.31) and (ii) we give a more detailed proof for existence of a unique 
" 7^ ''^-energy minimizer." 

In the following we denote by d{a, a') the Riemannian distance. The following formula 
is useful; 

d 

d{a,a') = sup{/(a) - f{a') \ f G C^{R',R), Y.^Vdf < l}- 

1=1 

Let a G R" be fixed and 0° = 0°(7) be as in (ESD for 7 G 

It is shown in [2] that, for any a', d{a,a') = inf{||7||-^ | 7 G } and infimum is 
actually attained. 

Lemma 9.2 If a and a' are sufficiently near, there exists 7 G K^' uniquely such that 
d{a,a') = W'jWh- (In other words, (A2) is satisfied.) 

Proof. In Lemma 4.2, p], it is shown that t ^ 0t(7) is a geodesies if and only if 
7t = vt for some v G R'^. For some open neighborhood f/ of G R'^(= TalV^), the map 
U 3 V ^ exp(v) defines a diffeomorphism onto its image, which is an open neighborhood 
of a'. Here, exp stands for the exponential map in the sense of Riemannian geometry. 
The path [0, 1] 3 t exp{tv) is the unique shortest geodesies from a to a' and is equal to 
t 0° = 0^(7) with 7t = vt. Consider a time change of this path: t 1— )■ 05j(t i)- Then, this 
is equal to t 0?(7) with 7^ = vR(t, 1) and d{a, a') = W'jWh (see the proof of Proposition 
4.3, [2]). 

It is sufficient to show that ||7||-h < for any other k G . First, we prove that, 

for any partition {0 = to < ti <■■■< ti} , it holds that 

I 

J2d{<Pl_^{k),<Pl{k))<\\k\\n. (9.1) 

i=l 

It is immediate from (19.11) that, unless t t— i- (f)f{k) is a time change of the geodesies 
1 1— )■ exp{tv), the inequality ||7||-h < \\k\\n holds. 
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For simplicity, we only show (19. ip when 1 = 2. The general case can be done in the 
same way. Let us recall some basic facts about Volterra kernel. For example, see Section 
2, [2], or Chapter 5, [15], and references therein. Set 

KH{s,t) = chs^'" [ {u- s)"~^u"-^du {t > s), 

J s 

where ch is a. certain positive constant. Its explicit value is not used in this paper. Denote 
by S the set of R'^- valued step functions on [0, 1]. Let 'H be the Hilbert space defined as 
the closure of S by the inner product 

{I[o^s]V,I[o,t]w)^ = R{s,t){v,w)jt,d, {t,s G [0, 1], v,w e R'^). 

The unitary isometry : ?^ — L^([0, 1], R*^) is given by 

{Ky){s) = J' f{t)^{s,t)dt. 

Since L'^{[0, 1],R°') is unitarily isometric to Cameron-Martin space T-L by an integral op- 
erator defined by 

{KHh){t)= [ KH{s,t)h{s)ds (/igL2([0,1],R'^)), (9.2) 

ii and "H are also unitarily isometric. 

Let < r < 1 and let k G Ti^ be such that 4>i{k) = a'. For simphcity, we will 
write zt = <ft{k). For any e > 0, there exist f,ge C^°°(R'^,R) such that Y^^iiViff and 
T.Uiy^9? < 1 and 

d{a, Zr) + d{zr, a') - 2e < f{zr) - f{a) + g{a) - g{zr) 

<E/ V,f{z,)dki + J2 V,g{z,)dki. (9.3) 
j=i j=i 

Assume k = Knh as in (I92D- Set Aj{s,t) = Vjf{zt)I{o<s<t<T} + Vjg{zt)I{r<t<ifi<s<t}- 
Then, by integration by parts, we have 

l\,f{Zs)dki = l\^f(z,)^l^^{s,t)hids)dt 

— — (s, t)Aj{s, t)hidsdt, 

{o<s<t<T} 

fv,g{zs)dkl = f\^g(zt)(^t^{s,t)hids)dt 

J T J T •-'0 

^ ^ '^s,t)Aj{s,t)hidsdt. 



{r<t<l,0<s<t} 9t 
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So, the right hand side of (19.31) is equal to 



E / ^-^{s,t)A,{s,t)^dsdt 

~[ J{0<s<t<l} 

= E [\f^{s,t)Ms,t)dt)hids 

j = l ^ 

<{tj2{f^('^t^\Ms,t)\dt)[f^^ 

< \\K'^l\\L2\\h\\L2 = \\l\\y^H\\h\\L2 = R{1, l)\\h\\L2 = ||/i||i2 = (9.4) 

Here, 1 stands for the constant function. Note that we used the unitary isometry of 
Hilbert spaces and the Schwarz inequahty 

d d I d I 

5]|A,(.,t)P,(.,OI< {$^|^,(^,t)P}'{ElA(^'^')r}'<l (forany.,t,tO. 

Thus, we have shown (19. ip for / = 2. The proof for general / is essentially the same. 

It is easy to see from (19. ip that, unless t zt = (pti^) is a time change of the geodesies 
t h-). exp(tw), it holds that d{a, a') < Let [0, 1] 9 t i— )■ G [0, 1] be a continuous and 

increasing function such that (Aq, Ai) = (0, 1) and A G Ti^i^lO, 1], R). Then, 9t := exp(Atf ) 
satisfies 9t = 0°(A.f ), or equivalently, 

d 

dOt = J2 VjiOtydXt, with do = a. 
i=i 

We have already seen that, for A = /?( ■ , 1), the equality d{a, a') = ||Af holds. It is 
sufficient to show that this equality fails for any other A. Such a A can be decomposed 
as A = i?( ■ , 1) + A, where A G '^■'^([0, 1], R) with Ai = 0. It is sufficient to show that 
i?( ■ , 1) and A are orthogonal in 'H-'^([0, 1], R), because it implies ||A||^h > ||i?( ■ , 1)||-^h. 
Let h G -^^^([0, 1], R) be such that Knh = A. From a well-known formula, 

R{t,T)= [ KH{t,s)KH{T,s)ds (t<T), 
Jo 
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we see that -R( ■ , 1) = Kh[Kh{1, ■ )]• Consequently, 

{R{-,l),X)nH = {Kh{1, ■),h)L2= [ KH{l,s)hsds = \i=0, (9.5) 

Jo 

which completes the proof of Lemma 19. 2[ I 

Next we will show (A3) when a and a' is sufficiently near. Combined with Lemma 
19. 2[ this lemma implies Proposition 19.11 

Lemma 9.3 Let 7 G K^' be the unique energy minimizer as in Lemma \9.S\ . If a and a' 
are sufficiently near, then (A3) is satisfied. (Note that a and a' may he nearer than in 
Lemma\EM) 

Proof. As in (the proof of) Lemma [9.2[ a' moves in a relatively compact neighborhood 
V contained in the geodesic coordinate centered at a. Recall that the correspondence 

V 3 a' ^ V ^ ^ = R{ - ,l)v e ^ iy. 

It is sufficient to show that, if V is taken small enough, then (A3)' is satisfied. Note that 
i? is uniquely obtained from 7 as in fl8.2l) - fl8.3l) . The only problem here is whether is 
small enough if V is small enough. 

First note that as a' varies in a relatively compact set, (pK'y) and Dcp^i^) = </>?(7)( ■ ) 
are bounded. Moreover, the positive constant c = 0(7) in (18.11) . with 7 replaced by 7, is 
bounded away from zero. 

By substituting D(j)^j{'y) (1 < j < d) for k, we obtain from (18.31) that 

((7, D<f>l,mH^, • • • , (7, D<PlMMC-' = {iy\...,iy'), 

where C stands for the covariance matrix in (18.11) . It is immediate from this representation 
that we can make \u\ sufficiently small by taking V is small enough. I 

Remark 9.4 Neither in nor in this appendix is a quantitative estimate of how near 
a and a' given explicitly. Hence, very precisely speaking, it is not completely shown that 
the result in [2] is included in ours. 

Remark 9.5 // we assume (Al) everywhere, then a Riemannian structure on R" is 
naturally induced as we explained above. If the case of the usual stochastic analysis (i.e., 
H = 1/2), (A2) and (A3) have a geometric meaning. (See Remark 3.2, 12^ . which was 
originally in J7^ \^.) First, ( A2) means that there is a unique shortest geodesies between 
a and a' . Second, (A3) or (A3)' means that these two points are not conjugate along the 
geodesies. So, Assumptions (Al)— (A3) are very mild and cover a lot of examples. 

It seems natural to guess from this that, in our case (i.e., 1/2 < H < 1), too. As- 
sumptions (Al)— (A3) are not bad. At this moment, however, the author does not have 
a nice example except the one in this Appendix. 



45 



References 



[1] Baudoin, F.; Hairer, M.; A version of Hormander's theorem for the fractional Brow- 
nian motion. Probab. Theory Related Fields 139 (2007), no. 3-4, 373-395. 

[2] Baudoin, F.; Ouyang, C; Small-time kernel expansion for solutions of stochastic dif- 
ferential equations driven by fractional Brownian motions. Stochastic Process. Appl. 
121 (2011), no. 4, 759-792. 

[3] Baudoin, F.; Ouyang, C; Gradient bounds for solutions of stochastic differential 
equations driven by fractional Brownian motions. Preprint. ' arXiv:1102.4601i 

[4] Bismut, J.-M.; Large deviations and the Malliavin calculus. Progress in Mathematics, 
45. Birkhauser Boston, Inc., Boston, MA, 1984. 

[5] Chronopoulou, A.; Tindel, S.; On inference for fractional differential equations. 
Preprint. ArXiv Math. 1104.3996. 

[6] Ciesielski, Z.; Kerkyacharian, G.; Roynette, B.; Quelques espaces fonctionnels as- 
socies a des processus gaussiens. Studia Math. 107 (1993), no. 2, 171-204. 

[7] Hu, Y.; Nualart, D.; Differential equations driven by Holder continuous functions of 
order greater than 1/2. Stochastic analysis and apphcations, 399-413, Abel Symp., 
2, Springer, Berlin, 2007. 

[8] Ikeda, N.; Watanabe, S.; Stochastic differential equations and diffusion processes. 
Second edition. North- Holland Publishing Co., Amsterdam; Kodansha, Ltd., Tokyo, 
1989. 

[9] Kuo, H.-H.; Gaussian measures in Banach spaces. Lecture Notes in Mathematics, 
Vol. 463. Springer- Verlag, Berlin-New York, 1975. 

[10] Lejay, A.; Controlled differential equations as Young integrals: a simple approach. J. 
Differential Equations 249 (2010), no. 8, 1777-1798. 

[11] Li, X.-D.; Lyons, T. J.; Smoothness of Ito maps and diffusion processes on path 
spaces. I. Ann. Sci. Ecole Norm. Sup. (4) 39 (2006), no. 4, 649-677. 

[12] Lyons, T. J.; Differential equations driven by rough signals. I. An extension of an 
inequality of L. C. Young. Math. Res. Lett. 1 (1994), no. 4, 451-464. 

[13] Lyons, T.; Qian, Z.; System control and rough paths. Oxford University Press, Ox- 
ford, 2002. 

[14] Molchanov, S. A.; Diffusion processes, and Riemannian geometry. Russian Math. 
Surveys 30 (1975), no. 1, 1-63. 



46 



[15] Nualart, D.; The Malliavin calculus and related topics. Second edition. Springer- 
Verlag, Berlin, 2006. 

[16] Nualart, D.; Ra§canu, A.; Differential equations driven by fractional Brownian mo- 
tion. Collect. Math. 53 (2002), no. 1, 55-81. 

[17] Nualart, D.; Saussereau, B.; Malliavin calculus for stochastic differential equations 
driven by a fractional Brownian motion. Stochastic Process. Appl. 119 (2009), no. 2, 
391-409. 

[18] Pratclli, M.; A remark on the variation of the fractional Brownian motion, 
Seminairc dc Probabilites Vol. 43, Lecture Notes Math. 2006 (2011). 215-219. 

[19] Shigekawa, 1.; Stochastic analysis. Translations of Mathematical Monographs, 224. 
American Mathematical Society, Providence, RI, 2004. 

[20] Watanabe, S.; Analysis of Wiener functionals (Malliavin calculus) and its applications 
to heat kernels. Ann. Probab. 15 (1987), no. 1, 1-39. 



47 



